Draft version March 8, 2013 

Preprint typeset using L^-T^X style cmulatcapj v. 5/2/11 



A SPECTROSCOPIC STUDY OF TYPE Ibc SUPERNOVA HOST GALAXIES 
FROM UNTARGETED SURVEYS 

N. E. Sanders 1 , A. M. Soderberg 1 , E. M. Levesque 2 , R. J. Foley 13 , R. Chornock 1 , D. Milisavljevic 1 , R. 
Margutti 1 , E. Berger 1 , M. R. Drout 1 , I. Czekala 1 , J. A. Dittmann 1 , 
Draft version March 8, 2013 

ABSTRACT 

We present the first spectroscopic study of the host environments of Type Ibc supernovae (SN Ibc) 
discovered exclusively by untargeted SN searches. Past studies of SN Ibc host environments have 
been biased towards high-mass, high-metallicity galaxies by focusing on SNe discovered in galaxy- 
targeted SN searches. Our new observations more than double the total number of spectroscopic stellar 
population age and metallicity measurements published for untargeted SN Ibc host environments. For 
the 12 SNe lb and 21 SNe Ic in our metallicity sample, we find median metallicities of 0.62 Zq and 
0.83 Zq, respectively, but determine that the discrepancy in the full distribution of metallicities is 
not statistically significant. This median difference would correspond to only a small difference in the 
mass loss via metal-line driven winds (< 30%), suggesting this does not play the dominant role in 
distinguishing SN lb and Ic progenitors. However, the median metallicity of the 7 broad-lined SN Ic 
(SN Ic-BL) in our sample is significantly lower, 0.45 Zq. The age of the young stellar population 
of SN Ic-BL host environments also seems to be lower than for SN lb and Ic, but our age sample 
is small. Combining all SN Ibc host environment spectroscopy from the literature to date does not 
reveal a significant difference in SN lb and Ic metallicities, but reinforces the significance of the lower 
metallicities for SN Ic-BL. This combined sample demonstrates that galaxy-targeted SN searches 
introduce a significant bias for studies seeking to infer the metallicity distribution of SN progenitors, 
and we identify and discuss other systematic effects that play smaller roles. We discuss the path 
forward for making progress on SN Ibc progenitor studies in the LSST era. 

This paper includes data gathered with the 6.5 m Magellan Telescopes located at Las Campanas 
Observatory, Chile. 

Subject headings: supernovae: general — galaxies: abundances — surveys 



1. INTRODUCTION 

Core-collapse supernovae show a diversity of absorp- 
tion features in their spectra near maximum light, re- 
flecting a diversity in the composition of the outer en- 
velope of their massive star progenitors at the end s 
of their lives (jFilippenkol 119971: IWooslev et all 120021) ■ 
In particular, some SNe show no evidence of hydro- 
gen (Type lb) or no evidence for either hydrogen or 
helium (Type Ic), suggesting extensive mass loss in 
the progenitor star sufficient to complete s tripping of 
the H and He layers of its outer envelope (lElias et al ' 
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1985t (Uomoto fe Kirshneri 119851: iClocchiatti et al l 119961 : 
Hachinger et al.l 120121 but see also lDessart et al.ll2012f ). 
Stellar evolutionary considerations point to two likely 
channels for these stripped-envelope core-collapse super- 
novae (Type Ibc supernovae) These channels are: (i) 
high-mass Wolf-Rayet (WR) stars with strong metal line- 
driven winds with rotation likely playing an important 
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role (|Wooslev et al.l [19951 iGeorgy et al.l [2011 . and (ii) 
lower-mass helium-stars in close-binary systems who lose 
their envelopes vi a Roche lobe overflow or c ommon en- 
velope evolutio n (Pods iadlowski et al.l 119921 lYoon et al.l 
[2010 : iDessart et al.ll2011| ). Searches for the progenitor 
stars of SN Ibc in pre-explosion imaging have yet to yield 
a progenitor detection, but have provided upper limits 
that challenge the hypothesis that their progenitors are 
massive WR stars like those seen in the Local Group 
(ISmarttll2009h . 

By measuring the metallicity of the host environments 
of Type Ibc SNe as a proxy for the metallicity of the 
progenitor stars, we may be able to distinguish between 
these two progenitor models. Mass loss in WR star s 
is enhanced at high metallicity (jVink fe de Koterll2005| ) . 
If the primary SN Ibc progenitor channel is single WR 
stars, then the rate of SN Ibc relative to SNe that show 
H features ( SN II) would be enhanced at high metallic- 
ity (see e.g. iPrantzos fc Boi ssicr 2003|) and the ratio of 
SN Ic to lb should similarly be higher. In binary progen- 
itor systems, massive primary stars may still strip their 
envelopes primarily via WR winds, but there exists an 
additional channel for SN Ic to be produced by relatively 
low mass (~ 12 Mq) stars which may dominate at low 
metallicity. This channel calls for the mass transfer to 
occur w hile the star is in t he core helium burning or later 
phases (jYoon et al.ll2010T ). Because either channel calls 
for massive, short-lived stars to produce the explosions, 
the metallicity of the SN host environment should be an 
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appropriate proxy for the mctallicity of the progenitor 
star. 

Additionally, a connection has emerged between long- 
duration Gamma Ray Bursts (long GRBs) and one par- 
ticu lar subtype of SN Ib c: b road-lined Type Ic SNe (Ic- 
BL: iKulkarni et al.l[l99l see lWooslev fc Blooml[2006l for 
a review). The broad and highly blueshifted absorption 
features of SNe Ic-BL indicate high phot ospheric expan- 
sion velocities, v p h > 2 x 10 4 km s _1 ()Iwamoto et al.l 
fl99l . This GRB-SN connection can be explained by 
the gravitational collapse of a massive (M > 20 Mq) 
progenitor star that produces a rapidly rotating and ac- 
creting compact object (central engine) that powers a rel- 
ativis tic outflow (the collapsar model, iMacFadven et al.l 
2001). However, radio observations demonstrate that 
only a small frac tion (< 1/3) of SNe Ic-BL ha rbor rela- 
tivists outflows (jSoderberg et al. 2006, 2010a). Because 
angular momentum loss due to metal line-driven winds 
could prevent the compact remnant from rotating fast 
enough to produce a relativistic outflow, a metallicity 
threshold (Z < 0.3 Z^) has been proposed for collapsars 
(jWooslev fc Hegerll2006h . 

Recent observations fundamentally challenge the rol e 
of metallicity in GRB production (jFrver et al.l 120071 ). 
SN 2009bb provides an example of a SN Ic-BL produced 
in a super-solar metallicity environment and harboring 
a cent ral engine (jLevesaue et alj|2010d: iSoderberg et al.l 
l2010a| ). and LGRBs with relatively high- metallicity host 
environments have now be en identified (jLevesque et al.1 
l2010bHGraham et al.ll2009h . In contrast, SN 2010ay was 
a SN Ic-BL with extreme explosion properties that oc- 
curred in a sub-solar metallicity environment and w ith- 
out evidence for a central engine (ISanders et al.ll201lT ; see 
also SN 2007bg. rYoTmg et alJl20ldj T~ 

In the past few years, several observational stud- 
ies have sought to measure t he character i stics o f the 
host environments of SN Ibc. iPrieto et al.l (|2008l) used 
SDSS spectroscopy to study the metallicity distribu- 
tion of 115 SNe (19 SN Ibc) and found that SN Ibc 
host environments are metal-enriche d compared to those 
of SN II. Extending the w o rk of iPrantzos fc Boissierl 
(l200l . 1 oissier &: Prantzosl ()2009[ ) reach a similar re- 
sult using SDSS photometry to esti mate metallici t ies for 
701 SN (98 SN Ibc) host galaxies. lArcavi et al.l pOlOh 
examine the host galaxies of core-collapse SNe discov- 
ered by the Palomar Transient Factory (PTF) and find 
that SN Ic are more common in high-metallicity en- 
vironments ("giant" host galaxies, M r < —18 mag), 
while SN lb, lib and Ic-BL dominate i n low-metallicity 
enviro nme nts ("dwarf" ho s t gala xies ) . lAnderson et al.l 
(|2010f ) and iLeloudas et al.l (|2011l) perform spectroscopy 
to measure the metallicity of the host environments of 28 
and 20 SN Ibc, respectively, finding no statistically sig- 
nificant difference between the mctallicit y distribution 
of SN l b and Ic. Ex t ending the work of IModiaz et al.l 
(|2008| ). IModiaz et al.l (|2011l ) performed a similar spec- 
troscopy study of 35 SN Ibc host environments, finding 
that SN Ic come from significantly higher-metallicity host 
environments tha n SN lb, with SN Ic -BL falling in be- 
tween Following iKellv et al.l ()2008l ). iKellv fc Kirshnerl 
(|2011f) have examined SDSS spectroscopy of 519 SNe (67 
SN Ibc), finding that SN Ic-BL preferentially occur in 
low-mass, low-metallicity host galaxies relative to other 
core-collapse SNe. 



However, there is a prominent observational bias that 
affects the progenitor metallicity distribution inferred 
from observations of most known SN host galaxies. The 
well-known relation between the luminosity and global- 
mctalli city of star-forming g alaxies (the L — Z relation, 
see e.g. iTremonti et al.ll2004l ) indicates that a metallicity 
distribution measured only from objects found in tar- 
geted SN searches, which look for transients in fields 
centered on nearby and luminous galaxies, will be biased 
towards high metallicities. 

Here we describe a new spectroscopic study of SN Ibc 
host galaxy metallicities unbiased with respect to the 
L — Z relation. We have obtained spectra of 60 host 
environments of SN Ibc discovered only by untargeted 
transient searches (13 SNe lb, 10 SNe lib, 24 SNe Ic, 8 
SNe Ic-BL, 3 of indeterminate type SN Ib/c, and 2 with 
AGN-dominated host environments, Section |2~3|) . Pre- 
vious spectroscopic studies of SN Ibc host environments 
have included relatively few SNe discovered by untar- 
geted searches, totaling < 40 objects. Our study doubles 
the existing sample of host environment spectroscopy 
for untargeted SN Ibc, offering considerable constrain- 
ing power for inferring the metallicity distribution of the 
parent population. 

In Section [3] we describe the characteristics of this 
sample, our optical observations, and our spectroscopic 
methodology. We present and analyze the host galaxy 
metallicities and other physical properties derived from 
these observations in Section [3] and combine with previ- 
ous spectroscopic surveys of SN Ibc host environments in 
Section 0J In Section [S] we discuss possible systematics 
affecting our results. In Section [5] we discuss our results 
in relation to SN Ibc progenitor models and the SN-GRB 
connection and we suggest implications for future studies 
of SN Ibc. We conclude in Section [71 

2. SAMPLE CONSTRUCTION 
2.1. SN sample 

We have observed the host galaxies of 60 SNe Ibc re- 
ported in the International Astronomical Union Circu- 
lars (IAUCsj^ and/or Astronomer's Telegrams^ between 
1990-2011. When transient searches operate by return- 
ing repeatedly to a prc-selectcd, typically bright, set of 
galaxies, we refer to the SNe found in those galaxies by 
those searches to be "targeted." We refer to any SN 
discovered by other means as "untargeted," including 
discoveries by wide-field optical surveys, SNe identified 
by targeted searches in anonymous background galaxies, 
and SNe discovered serendipitously during observations 
of unrelated objects. We observed only untargeted dis- 
coveries and prioritized those SNe with reliable classifica- 
tions, host galaxies which did not already have previously 
published metallicity measurements, and which were vis- 
ible at low airmass during the time of our observations. 

In total, we present optical spectroscopy with S/N suf- 
ficient for metallicity measurements for 50 host galax- 
ies, with median redshift z = 0.036. The discoverer of 
each SN in our sample is listed in Table [T] Spectro- 
scopic metallicity estimates have been previously pub- 
lished for only 15 of these galaxies. A comparison 

5 fl.html 

6 http:/ /www. astronomerstelegram.org/ 
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to previous spectroscopic studies of SN Ibc host envi- 
ronments is presented in Section 14.11 additionally, the 
host environment of SN 2009j f was previously s tudied by 
IValenti et afl (|2Qlll ) (see also lSahu et~aT1l2011| ). 

For those SNe in our sample whose classification spec- 
troscopy is not well discussed in the circulars or the lit- 
erature, we confirm the spectroscopic classification by 
contacting the authors of the classification circular to re- 
evaluate their origi nal spectra using e.g., th e Supernova 
Identification tool (|Blondin fe T onrv 2007|) (see private 
communications in Table [1]). We construct a "Gold" 
sample of 48 SNe for which the classification spectrum 
suggests a clear subtype (but see Section I5.3|) ; a "Sil- 
ver" sample of 10 SN where the classification is based 
on a spectrum with poor S/N. SNe where we are not 
able to distinguish between two sub-types are listed as 
"Undetermined Ibc." For the purposes of this study, we 
consider the peculiar objec t SN 2011hw (classified as 
Type Ibn; iSmith et al.ll2012l ) to be of indeterminate type, 
we consid er SN 20061c (whose He lines were weak) a 
Type lb (jLeloudas et al.l |2011[ ). and we group the pe- 
culiar Type lib SNe 2002gz and 2010cn together with 
the other Type lib SNe. 

2.2. Spectroscopic Observations 
We obtained acquisition images in r-band and long-slit 
spectra (~ 3600 - 8500 A) of 25 SN Ibc host galaxies us- 
ing the Low Dispersion Survey Spectrograph 3 (LDSS3) 
instrument on the Magellan-Clay Telescope, 23 using 
the BlueChannel Spectrograph (BC) on the MMT (for 
which there is no associated imager) , 9 using the Inamori- 
Magellan Areal Camera & Spectrograph (IMACS) in- 
strument on the Magellan-Baade Telescope, and 1 using 
the Gemini Multi-Object Spectrograph (GMOS) instru- 
ment on the Gemini-North Telescope. The host galax- 
ies were observed at the parallactic angle, except for 
some IMACS spectra which were obtained with the at- 
mospheric dispersion compensator. 

When possible, the spectrum was extracted at the loca- 
tion of the SN explosion site within the host galaxy. We 
consider the spectrum to sample the explosion site ( "E" ; 
43 SNe) when we extract at the position of the explo- 
sion site and the slit width corresponds to a physical size 
< 2 kpc (z < 0.11 for a 1" slit), and otherwise consider 
the spectrum nuclear ("N"; 15 SNe, see Table EJ. Explo- 
sion site spectroscopy represents a luminosity-weighted 
average of the physical properties of the star-forming gas 
within the extracted region. However, because our SN 
sample consists of events discovered by galaxy-impartial 
surveys, the host galaxies are typically smaller and more 
distant than those observed in previous studies. In some 
cases, the explosion sites of SNe in the outskirts of their 
host galaxies did not have sufficient flux for spectroscopy 
and instead we extract at the galaxy nucleus. Some 
intrinsically-dim host galaxies have apparent sizes so 
small that a significant fraction of all the galaxy light 
will fall in the slit, even at relatively-low redshift (see Fig- 
ur e [H . For compariso n, the explosion site spectroscopy 
of iModiaz et al.l (120111 ) includes host galaxies at a max- 
imum redshift of z = 0.183 using a 1" slit (~ 3 kpc) and 
the st udies of lPrieto etall (|2008h and lKellv fc Kirshnerl 
(120111 ) employ 3" SDSS fiber spectroscopy for SN Ibc at 
z < 0.04 (~ 2 kpc) and 0.07 (~ 4 kpc), respectively. 



2.3. Spectroscopic analysis 

We employed standard two-dimensional long-slit im- 
age reduction and spectral extraction routines in IRArfl 
Examples of reduced spectra are displayed in Figure [2J 

The line fluxes of prominent nebular emission lines (Ha 
and H/3; [O II] A37270, [O III] A4363, [O III] A4959, 
[O III] A5007, [N II] AA6548, 6584; and [S II] AA6717, 
6731) were measured by fitting Gaussian functions 
to their profiles usi ng a Markov Ch ain Monte Carlo 
(MCMC) technique (|Patil et all 120101 ) . We fit the pro- 
files to a wavelength range 20 A in width centered on 
the rest frame wavelength of each line. We fit a linear 
continuum to 20 A regions of the spectra off the wings 
of each line and use this continuum measurement to cal- 
culate equivalent widths for spectral lines. We constrain 
the amplitude of the Gaussian to be positive and fit a 
single redshift and line width for the set of Balmer lines 
and an independent redshift and width for the forbidden 
lines. We adopt the fitted redshift of the Balmer lines 
as the redshift of the host galaxy. We obtain estimates 
of the uncertainty in the line flux from the MCMC trace 
and require a detection confidence of 999c0. 

The line fluxes measured for each host environment are 
presented in Table [3l We use observations of standard 
stars from the same night as the host environment obser- 
vations to achieve relative flux calibration. No correction 
for underlying stellar absorption has been made, which 
could potentially affect the measurement of the Balmer 
line fluxes (particularly H/3) for galaxies with significant 
stellar continuum flux. Only 5 of 58 spectra (those of 
SN 2002ex, 2006ip, 20071, 2007A, and PTFlOvgv) have 
stellar continuum levels that may indicate significant un- 
derlying absorption, and we therefore do not measure 
quantities that depend on H/3 for these objects, includ- 
ing extinction (Ay). 

We tested f or A GN contamination using a 
I Baldwin et all (|1981[ ) excitation-mechanism diag- 
nostic diagram and the classification scheme of 
(jKauffmann et all 12001 . We found emission lines 
from the host galaxies to be consistent with typical 
star-forming galaxies with few exceptions. The host 
environments of SN 1991D and 2003jp were excluded 
from our sample (and all analysis below) due to sig- 
nificant AGN contamin ation that would bias certain 
metallicity diagnostics (jKewlev fc E llison 20081) . We 
neglect the potential effects of AGN contamination 
for SN 2006jo (Type lb) and PTFlOxla (Ib/c), whose 
host environments show evidence for a composite 
classification. Some objects have supernova flux con- 
tamination (2010ay, 2011bv, 2011gh, 2011ip, LSQ11JW, 
PTFlOaav), although the broad SN absorption features 
do not affect the measurement of the flux in the narrow 
galaxy emission lines. 

7 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 

8 The [O II] AA3726, 3729 doublet is not resolved in these spectra, 
and we refer to the sum we effectively measure as [O II] A3727. 

9 Assuming a normal distribution, the 99% confidence interval 
corresponds to the requirement that the median value be 2.576cr 
greater than zero. Here we calculate a as the difference between 
the 50th and 16th percentile values of the amplitude distribution 
in the MCMC trace. 
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TABLE 1 

Properties of SNe in Untargeted Sample 



SN 


Type a 


Samplc b 


z 


Slit width (kpc) c 


Discoverer d 


1991R 


Ibc 


S 


0.035 


0.7 


TVTAT 1 j C TV f T/ 

McN aught fe McKenzic 


2002cx 


lb 


G 


0.037 


1.1 




zuuzgz 


lib-pec 


c 
5 


U.Uou 


i i 




ZUUoCV 


IC 


C 7 


n no/i 
U.UZ4 


U.o 




ZUU4CI 


ID 


C 

O 


U.Z4o 




rucllo ct al. (zUU4b) 


ZUU41D 


IC 




U.Uob 


1 1 
1 . 1 


oiJoo 


zUUonm 


lb 


It 


n riQ/i 
U.Uo4 


1 n 
l.U 


oJJoo 


zUUonb 


lc-r3L 


Ct 


0.024 


0.5 


ymrnby ct (^Uub; 


zUUbip 


ic 


Ct 


0.031 


0.6 


^^^^^^^ 


200bir 


Ic 


G 


0.021 


0.4 


SNF 


200bJO 


lb 


Ct 


0.077 


1.4 


nrvnn 
bJJbb 


ZOOblC 


ID 


It 


0.016 


0.3 


bJJbb 


2006nx 


Ic-BL 


G 


0.137 


2.4 


SDSS 


2006tq 


Ic 


S 


0.261 


3.0 


ESSENCE 


20071 


Ic-BL 


G 


0.022 


0.4 


LOSS 


2007az 


lb 


G 


0.035 


1.0 


LOSS 


2007br 


lib 


G 


0.053 


1.0 


SNF 


2007cc 


Ic-BL 


S 


0.046 


0.9 


Ouimbv et al. (20071 


2007db 


Ic 


G 


0.048 


0.7 


SNF 


2007oa 


lib 


G 


0.040 


0.8 


SNF 


2007fT 


Ic 


S 


0.049 


0.7 


SNF 


2007gg 


lb 


G 


0.038 


0.8 


SNF 


2007gl 


Ic 


G 


0.028 


0. 1 


SNF 


2007hb 


Ic 


G 


0.021 


0. i 


SNF 


2007hl 


Ic 


S 


0.056 


1.1 


SNF 


2007hn 


Ic 


G 


0.028 


0.6 


SNF 


2008ao 


Ic 


G 


0.015 


0.3 


Dimai ("2008") 


2008fi 


lib 


G 


0.026 


0.5 


Skvarc & Mikuz (2008) 


2008gc 


lb 


G 


0.049 


0.7 


CHASE 


2008ik 


Ic 


G 


0.064 


0.9 


CHASE 


2008im 


lb 


G 


0.008 


0.2 


Oksanen (2008) 


2008iu 


Ic-BL 


S 


0.130 


1.6 


CRTS 


2009hu 


lb 


G 


0.117 


2.1 


Sand et al. (2009") 


2009jf 


lb 


G 


0.008 


0.1 


LOSS/PTF 


2009nl 


Ic 


Ct 


0.113 


1.8 


CRTS 


2010Q 


Ic 


G 


0.054 


1.0 


CRTS 


2010ah 


Ic-BL 


G 


0.050 


1.0 


PTF 


2010am 


lib 


G 


0.020 


0.4 


CRTS 


2010ay 


Ic-BL 


G 


0.067 


1.3 


CRTS 


2010cn 


Ib/IIb-pec 


G 


0.026 


0.5 


CHASE 


20101z 


Ic 


G 


0.090 


1.2 


CRTS 


2011D 


lib 


G 


0.023 


0.3 


CRTS 


2011V 


lib 


G 


0.014 


0.3 


CRTS 


2011bv 


lib 


G 


0.072 


1.0 


CRTS 


2011cs 


Ic 


G 


0.101 


1.8 


CRTS 


2011gh 


Ib/c 


S 


0.018 


0.4 


CRTS 


2011hw 


Ibn 


G 


0.021 


0.4 


Dmtmiana ct al. (2011J 


2011ip 


Ic 


S 


0.051 


1.0 


Denisenko et al. f2011 , 


2011it 


Ic 


G 


0.016 


0.3 


Ciabattari et al. (2011) 


LSQ11JW 


lb 


G 


0.020 


0.4 


LSQ 


PTF09dfk 


lb 


G 


0.016 


0.2 


PTF 


PTF09dxv 


lib 


G 


0.032 


0.4 


PTF 


PTF09iqd 


Ic 


G 


0.034 


0.7 


PTF 


PTF09q 


Ic 


G 


0.090 


1.7 


PTF 


PTFlOaavz 


Ic-BL 


G 


0.063 


0.9 


PTF 


PTFlObip 


Ic 


G 


0.051 


1.0 


PTF 


PTFlOvgv 


Ic 


G 


0.015 


0.3 


PTF 


PTFllhyg 


Ic 


G 


0.028 


0.6 


PTF 



Classification 



Jje^bjindgm^et^L 
S. 

Hamuv et al. (2003, 
^^^^^^^^^_CJiornock 

M. T. Botticella 
G. Leloudas 
Leloudas et al. (2011) 
Modiaz etld"72011) 
Modiaz et al. (20TT) 
Leloudas et al. (2011) 
Leloudas et al. (2011) 
Blondin et al. (2006) 
Modiaz et al. f20PH 
Challis et al. (2007) f 
Modiaz et al. (2011) 
Silverman et al. (2007) 



R. Thomas 
Odewahn & 



Tnomas 
Thomas 
Thomas 
Thomas 
Thomas 
Thomas 
Thomas 



Terrazas (2007 



Leloudas et al 




A. Drake 
A. Drake 
Corsi et al. (2011) 
Ccn ko et al. f^oTo) 
Sand ers et al. (2011) 



A. Drake 

Marion et al. (2011) 

Milisavlievic et al. (2011) 

Prieto (20lTl 

Drake et al. (2011b) 

Magill et al. (2011) 

Smith et al. (2012) 

S., S. Valenti 

S. 

S. 

A. Gal- Yam 
A. Gal- Yam 
A. Gal- Yam 
A. Gal- Yam 
S. 

A. Gal- Yam 
Corsi et al. (2012) 
A. Gal- Yam 



a A detailed discussion of SN classification is given in Section 12. II 

k G indicates the Gold sample and S indicates the Silver sample, as defined according to security of spectroscopic classification in Section l2. II 
c The size of the spectroscopic slit (sec Table [2l in physical units at the distance of the SN. 

Reference for the discovery of the supernova. Acronyms for untargeted SN searches arc as follows: Catalina Real-time Transient Survey 
(CRTS; IDrake et al~ll2009a!) , the Equat ion of State: SupErNo vae trace Cosmic Expansion program (E SSENCE; IMiknaitis et aLlfeOOTT) , the La 
Silla-QUEST Variability Survey (LS Q; [H adjiyska c t a.1.1120 1 lT> ■ the Palomar Transient Factory (PTF; ILaw et al1l2009l) . the Sloan Digital Sky 
Survey-II Supernova Survey (SDSS; IFricman ct al.H2008l) . and the Nearby S upcrnova Factory (SNF, including supcrnovae discovered by NEAT; 
Aldcring ct al. 2002). Some objects were discovered serendipitously in background galaxies during targeted SN searches or other wide-field 
surveys, including the CHilean Automatic Supernova sEarch (CHASE; IPignata et al.H2009H and the Lick Observatory Supernova Search (LOSS; 
ILi et al.|[20lll) . 

c Reference for spectroscopic classification of SN type. Author names refer to private communications with observers who performed spectroscopy 
of the SN. "S." indicates our own spectroscopy. 

^ Additionally, private communication with S. Blondin (2012) indicates that the classification spectrum for SN 2006tq is not of S/N sufficient 
for the Gold sample. 

s Varying spectroscopic classifications of SN 2010cn arc reported by[Marion et all (f20To|) ; |Silverman et al.| i2010fi ;|Folalclli & Stritzingcr (20HJ); 
we consider it a SN lib for the purposes on our anaalysis. 
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TABLE 2 

SN Host Galaxy Sample and Observing Configurations 





SN Ck2000 


SN 62000 


Date (UT) 


Instrument 


Disperser a 


Exp. time (s) b 


Slit width 


E/N c 


1991D^ 


1 Q./l 1 .1 Q C.Q. 
lO.'il.lO.OO 


-l^.OO.1 1 .u 


ZUUo 


iviay ou 


LiJOOO 


VPH_A11 


UUU 


10" 
± . u 


\? 


1991R 




+ 1 9 00-43 9 


ZUUo 


iviay ou 


T.DSS3 


VPH.A11 


600 


1.0" 




ZUUZcX 


99-0.0-0.0 7Q 

zz. uy.uu. . y 


in.Qfi.Ot; Q 
1U.OU.ZO.O 


onnc 

ZUUo 


T„ nQ on 

j unc zy 


LiJODO 


V 1 1 1 1 ' 1 lie / rveo 


1 900 /I 090 
1ZUU/ J-UZU 


1 e \" 

1.0 


\? 


90090-7 
zuuzgz 


UZ.O'l.lU.OU 


nn. c;'3. 1 a 

UU.OO. IO.Z 


onns 

ZUUo 


j unc zy 


TA/r Ar^ 


Qnn 17 t; 


1 900 
1ZUU 


7" 
u. 1 




ZUUOcV 


19.1 n.QI on 
XO.XU.oX.C5U 


01 -^Q-zlQ R 

— zi ,oy .4y .d 


onno 
ZUUo 


June 1 


IjlJooO 


VPLT All 


QOO 

yuu 


10" 

1 . u 


L' 


2003jp 


23:26:03.28 


-08:59:22.7 


2007 


Dec. 17 


LDSS3 


T rT^TT All 

VBll.All 


600 


0.75 " 


N 


2004cf 


14:11:05.77 


-11:44:09.4 


2006 


June 30 


LDSS3 


VPH_Rcd 


1800 


1.5" 


N 


2004ib 


02:40:56.40 


r\r\ 10 Ad 

-00:10:48.3 


2007 


Dec. 16 


LDSS3 


t rriTT An 


1200 


1.0" 


E 


2005hm 


21:39:00.65 


O 1 o 1 OO 'I 

-01:01:38.7 


2006 


June 30 


LDSS3 


V BH.Bluc / Hcd 


01 c\r\ /1 f OA 

2100/1500 


1.5" 


E 


2005nb 


12:13:37.61 


+ 16:07:16.2 


2008 


May 31 


LDSS3 


T 7"T»TT All 

V-PH_A11 


415 


1.0" 


E 


2006ip 


OO .11) Ml . ' _ } 

23:48:31.68 


00 00 p" t 

-02:08:57.3 


2007 


Dec. 16 


LDSS3 


t rr>TT An 


1200 


1.0" 


N 


2006ir 


23:04:35.68 


+07:36:21.5 


2007 


Dec. 14 


LDSS3 


t rr>TT An 


1200 


1.0" 


E 


2006jo 


/"\ 1 1.1 TO 

01:23:14.72 


-00:19:46.7 


2006 


Dec. 24 


LDSS3 


VPH_Bmc/Rcd 


1000/600 


1.0" 


N 


20061c 


22:44:24.48 


-00:09:53.5 


2007 


Dec. 14 


LDSS3 


VPH.A11 


900 


1.0" 


N 


2006nx 


03:33:30.63 


-00:40:38.2 


2006 


Dec. 24 


LDSS3 


VPH_Bluc / Red 


1500/750 


1.0" 


N 


2006tq 


02:10:00.70 


+04:06:00.9 


2012 


Jan. 18 


LDSS3 


VPH.ALL 


1200 


0.75" 


N 


20071 


11:59:13.15 


-01:36:18.9 


2008 


Jan. 17 


LDSS3 


~\ T T~\ T T All 

VPH_A11 


1800 


1.0" 


N 


2007az 


08:25:23.80 


+69:54:29.4 


2011 


Dec. 24 


BC 


300GPM 


1200 


1.5" 


E 


2007br 


11:15:39.93 


-04:22:47.8 


2007 


Dec. 14 


LDSS3 


VPH.A11 


1200 


1.0" 


N 


2007cc 


12:10:17.96 


+48:43:31.5 


2012 


Jan. 01 


BC 


300GPM 


1800 


1.0" 


E 


2007db 


11:17:10.30 


-06:11:48.6 


2007 


Dec. 15 


LDSS3 


t rnTT All 

VPH_A11 


900 


0.75 " 


N 


2007ca 


IP" fO /I /"> O *~7 

15:53:46.27 


^7 oo 1 r- f 

-27:02:15.5 


2008 


May 30 


LDSS3 


t rr>TT An 

VPH_A11 


1200 


1.0" 


E 


2007n 


01 o .4 i o ri i 

01:24:10.24 


1 00 00 a r\ p" 

+09:00:40.5 


2007 


Dec. 17 


LDSS3 


t rriTT An 

VPH_A11 


1200 


0.75" 


N 


2007gg 


00:28:12.51 


+00:07:04.8 


2011 


Sep. 06 


BC 


300GPM 


1200 


1.0" 


E 


2007gl 


01 o c t; r\ 11 

01:05:50.11 


1 r\r\ r\o a 1 


2011 


Sep, 21 


T A K A 

IMALo 


oOU-17.5 


1200 


0.7" 


E 


zUU t riD 


UO.Uz.Ul.zo 


01 .n?. 1 
— zl.UY .00.1 


2011 


Inov. ly 


TA/T A f^C 
IIVIAL^O 


Qnn 17 c; 
oUU-1 ( .0 


ftnn 
DUU 


no" 
u.y 


\r 

IL 


ZUU i ill 


on-c,n-n7 ic\ 
zu.OU.U 1 . i U 


— Ul.Oo.ou.4 


2008 


May 30 


T F1QC3Q 
LlJooo 


A/PTT All 


onn 

yuu 


1 n " 
l.U 


c 1 
Ji 


ZUU / illl 


Zl.UZ.4D.oO 


n/i ■n^.of; 

— U4.UD.ZD.Z 


2008 


May 30 


T FlQ^Q 
IjIJooo 


A/PTT All 
V r H_rVll 


DUU 


1 n " 
l.U 


IN 


zUUoao 


Uo.U ( .4D.DD 


1 QQ -OO-nfi 


2012 


Jan. 01 


jL>1 


QnnPDi\j 
oUULtI^IVI 


1 Qnn 
loUU 


1 n " 
l.U 


\r 


zuuon 


m - c,q.oq 1 7 
Ul.Oo.zo.l f 


_i_on-oi .oo a 
-|-zy .zi .zo.4 


2012 


Jan. 19 




Qnnr^T3i\/r 

OUULrl^lVl 


1 onn 
IzUU 


1 n " 
l.U 




zuuogc 


no-i n-QK fiQ 
Uz.lU.oD.Do 


— 0o.40.0y .0 


2012 


Jan. 18 


IjIJooo 


A/PTT All 
V r H_rVll 


finn 
DUU 


U. 1 




ZUUolK 


nQ-Qfi-no ka 
Uo.oD.uy.04 


qc.i q. nn 7 
— oO. lo.UU. i 


om 
zUlz 


Jan. zu 


LlJooo 


upu All 
V r rl_All 


finn 
DUU 


n 7c; " 
U. ( 


\r 

IL 


ZUUollTl 


ha .m -no 1 c. 
U4.U1 -Uz. 10 


-p ( 4.U0.4O.0 


2011 


Dec. 24 




OUULrl^lVl 


ouu 


1 K " 
1.0 


c 1 
Ji 


ZUUolU 


H/1 -Qfi. c,c on 
U4.oD.D0.zU 


nn-oi .Qt; k 
— UU.zl.oO.D 


2011 


Sep. 21 


llVl/iOO 


Qnn 17 c, 
oUU-1 1 .0 


1 onn 
IzUU 


n 7" 




zUUynu 


1 /i . cc q . on QO 
14.0o.zy.oz 


1 1Q.QC.Q1 1 

+ lo.o0.ol.l 


2012 


Jan. 01 


"DP 


OUULrr 1V1 


1 onn 
loUU 


1 n " 
l.U 


1 

iL 


zuuyji 


OQ.n/i -co ciq 
zo. U4.0z.yo 


1 1 o. 1 n. en c 
-f-lz.ly.Oy.O 


2011 


Nov. 18 


IIVIAL^O 


Qnn 1 7 e 
oUU-1 1 .0 


enn 
DUU 


n 7" 


li 


zuuym 


Uo.oy.4 1 . / 


1 1 ■ 1 q ■ k n 

—11 . lo.zo.U 


om 1 
ZU11 


oep. ZU 


TA/T A 
llVl/iOO 


Qnn 17 c, 
oUU-1 1 .0 


1 onn 
IzUU 


no" 
u.y 


il^ 


om nn 
zUlU*^ 


1 n-OR-o'v 1 1 
lU.zD.z ( .11 


1 Qn.m ■ cn n 
H-oy.U1.0U.y 


2012 


Jan. 18 


jL>1 


oUUt-j-r 1V1 


1 Qnn 
loUU 


1 n " 
l.U 


\r 

IL 


2010ah 


11:44:02.99 


+55:41:27.6 


T 




BC 


300GPM 


1 t" r\r\ i 1 ovir\ 

1500+1200 


1.0" 


E 


2010am 


09:33:01.75 


+ 15:49:08.8 


2012 


Jan. 01 


BC 


300GPM 


1800 


1.0" 


E 


2010ay 


12:35:27.19 


+27:04:02.8 


2010 


Apr. 11 


GMOS 


R400 


1800 


1.0" 


E 


2010cn 


11:04:06.57 


+04:49:58.7 


2012 


Jan. 18 


BC 


300GPM 


1200 


1.0" 


E 


20101z 


01:50:20.32 


-21:44:31.9 


2012 


Jan. 19 


LDSS3 


i r"nTT A11 

VPH_A11 


1500 


0.75" 


E 


2011D 


03:02:14.58 


+17:20:58.9 


2011 


Nov. 29 


IMACS 


300-17.5 


500 


0.7" 


E 


2011V 


09:27:38.76 


+28:47:27.2 


om 
zUlz 


Jan. 01 


BC 


300GPM 


1800 


1.0" 


E 


201 lbv 


13:02:53.57 


-04:02:36.0 


2012 


Jan. 19 


LDSS3 


T mTT All 

VPH_A11 


1200 


0.75" 


E 


2011cs 


12:08:01.08 


-f-^y. 10.00. u 


2012 


Jan. 19 


DP 


Qnnr^PA/T 
OUUVjlr^lVl 


1 ^00 

10UU 


10" 

l.U 


u- 


201 lgh 


03:16:54.20 


+25:54:14.6 


2012 


Jan. 18 


BC 


300GPM 


600 


1.0" 


E 


201 lhw 


22:26:14.54 


+34:12:59.1 


2011 


Nov. 29 


BC 


300GPM 


600 


1.0" 


E 


2011ip 


1:13:47.59 


-12:41:06.0 


2011 


Dec. 31 


BC 


300GPM 


1800 


1.0" 


E 


201 lit 


22:02:44.45 


+31:41:49.1 


2011 


Dec. 31 


BC 


300GPM 


600 


1.0" 


E 


LSQ11JW 


02:04:47.40 


+00:50:06.0 


2011 


Dec. 31 


BC 


300GPM 


1200 


1.0" 


E 


PTF09q 


12:24:50.11 


+08:25:58.8 


2012 


Jan. 18 


BC 


300GPM 


600 


1.0" 


E 


PTF09dfk 


23:09:13.42 


+07:48:15.4 


2011 


Sep. 21 


IMACS 


300-17.5 


900 


0.7" 


E 


PTF09dxv 


23:08:34.73 


+18:56:13.7 


2011 


Nov. 18 


IMACS 


300-17.5 


500 


0.7" 


N 


PTF09iqd 


02:35:23.23 


+40:17:08.7 


2012 


Jan. 18 


BC 


300GPM 


500 


1.0" 


N 


PTFlOaavz 11:20:13.36 


+03:44:45.2 


2011 


Jan. 13 


LDSS3 


VPH.A11 


1200 


0.75" 


E 


PTFlObip 


12:34:10.52 


+08:21:48.5 


2012 


Jan. 01 


BC 


300GPM 


1800 


1.0" 


E 


PTFlOvgv 


22:16:01.17 


+40:52:03.3 


2011 


Sep. 06 


BC 


300GPM 


1200 


1.0" 


E 


PTFllhyg 


23:27:57.34 


+08:46:38.0 


2011 


Sep. 06 


BC 


300GPM 


1200 


1.0" 


E 



a The grism or grating used in the spectroscopic observation. When multiple configurations were used, we list both, separated with 
a slash. 

The exposure time of the spectroscopic observations. When two times are given separated by a slash, they correspond to two 
different dispcrscrs. 

c E indicates a slit position on the SN explosion site, N indicates slit placement at the galaxy center (sec Scction |2.2t 
^ Objects excluded from the sample due to AGN contamination; sec Scction |2.3l 
^ Spectra acquired on 2011 May 06 and 2012 Jan. 18 were coadded. 



6 



Sanders et al. 



a) NGC 7479 / SN 2009jf 
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20.O" 


-b) PTFIOaavz 
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Fig. 1. — Example of optical spectroscopy procedure for two SN Ibc host galaxies in our sample. Shown are the r-band acquisition images 
(top) and 2D-spectrum (bottom) obtained using Magellan/IM ACS and LDSS3 (see Table[2j| for a) the Type lb SN 2009jf host environment 
in the galaxy NGC 7479 (z = 0.008) and b) the Type Ic-BL PTFIOaavz (z = 0.063). The slit (blue rectangle) was aligned with the SN 
explosion site (red circle). The IMACS chip gap is visible in the spectrum at ~ 4100 A. 
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c) 
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2007ce 

Ic-BL z=0.046 ii i 
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d) 



2008fi 
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Rest wavelength (A) 

Fig. 2. — Examples of spectra from our SN Ibc host environment observations. The panels are labeled with the names, classifications, 
and redshifts of the SNe and represent a) a high-S/N spectrum dominated by nebular emission from the host environment, b) a spectrum 
dominated by SN flux, c) a spectrum dominated by stellar continuum, and d) a low S/N spectrum dominated by nebular emission from 
the host environment. The gray vertical lines illustrate the position of the major nebular emission features used in the metallicity analysis. 
The continuum level can be judged against the dashed horizontal line, illustrating zero flux. 
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3. HOST ENVIRONMENT PROPERTIES 

From our optical spectra we measure the magnitude 
of dust extinction, metallicity, young stellar population 
ages, and Wolf-Rayet star populations of the host envi- 
ronments of the SN Ibc in our sample. We discuss each 
of these in the following sections. 

3.1. Dust extinction 

We estimate the line-of-sight extinction for each SN 
host galaxy spectrum (including both Galactic and in- 
trinsic reddening) from the Balmer flux decrement. We 
assume Ha/H/3 = 2.85, which corresponds to T = 
10,0 00 K and n K = 10 4 cm" 3 f or Case B recombina- 
ti on (lOsterbrock fe Ferlandl 120061 ). The extinction curve 
of lCardelli et al.l (|1989l ) was applied to correct individual 
line fluxes for reddening, assuming Ry =3.1. The value 
of the extinction derived for each host galaxy is presented 
in Table El 

In Figurc|3l we show the distribution of visual-band ex- 
tinction (Ay) derived for the SN host environments from 
the Balmer decrement, for all spectra (solid lines) and for 
only explosion-site spectroscopy (dashed lines). Galactic 
extinct ion has been subtracted usin g the infrared dust 
map of iSchlaflv fc Finkbeinerl (|2011l ). The difference in- 
troduced by restricting the sample to explosion sites is 
small. 

Using the Kolmogorov-Smirnov (KS) test, we do not 
identify a statistical difference between the extinction 
distributions of SN lb and Ic (with KS p-value, j»ks = 
0.85). Combining the 12 SN lb and 23 SN Ic with A v 
measurements and subtracting Galactic extinction, the 
median and 14th and 86th percentile values (lcr) arc 
A v = 0.2 

-o% ma g- The median extinction values for 
the SN lib and Ic-BL are consistent with 0, but given 
the small sample sizes (N ~ 9 and 6, respectively), the 
difference from the SN Ib+Ic distribution is not signifi- 
cant (p KS > 0.1). 

In Figu re [3] we also show the A y distribution mea- 
sured by IKellv fc Kirshnerl (pOll for SN II and for 
the combination of SN lb and Ic. A consistent ex- 
tinction dist r ibutio n was estimated for 19 SN Ibc by 
IDrout et al.l (1201 ID based on l ight curve colors. For 
SN Ibc. IKellv fc Kirshnerl (pOlTI ) find median and lcr val- 
ues of Ay = I.IIq 6 mag, with smaller values for SN II. 
We find Ay values with a more extreme range, with 
some SNe having no measurable redenning and some as 
large as Ay pa 5 mag. This discrepancy is likely due 
to methodol ogical differences. Becau se they employ 3" 
SDSS fibers, IKellv fc Kirshnerl ((2011 probe gas in the 
nuclear region of the host galaxy, whereas most of our 
spectra are from the SN explosion site where the line of 
sight could be significantly different. However, in some 
cases the discovery magnitudes reported for the SNe do 
not allow for extinctions as large as we measure for the 
host environment- This implies that the line of sight 
to the SN may be substantially different than that repre- 
sented by the integrated light from the star forming gas 
we observe. 

]°e.g. SN 2011D, discovered at m = 18.2 mag by INarla etTaD 

1201 U . implying M < —20 adopting the Balmer decrement extinc- 
tion 




A v (mag) 



Fig. 3. — Cumulative distribution of visual-band extinction (Ay) 
measured for host environments of SN of different types from the 
spectra collected in our sample. The dashed lines illustrate the 
same distributions, for only explosion-site spectroscopy. The num- 
ber of objects included in each sample is given in the legend, with 
explosio n-site only number s in pa rentheses. The distributions mea- 
sured by Kcllv & Kirshncr (2011) for SDSS fiber spectra are shown 
in black. 

3.2. Metallicity Estimation 

Host galaxy oxygen abundances were calculated from 
the extinction-corrected line flux ratios according to a 
number of independently calibrated abundance diagnos- 
tics, described below. We report the abundance derived 
from each method in Tabic 01 

We derive "direct" oxygen abundance estimates by 
estimating the electron temperature of the gas' dom- 
inant excitation zone, which is only possible if the 
O III] A4363 line is detected. Following iLevesque et al.1 
2010a| ). the electron temperatures are estimated us- 
ing IRAF's five-level n ebular modeling package nebular 
(|Shaw fc Dufourl [19941) . The nebular task temden is 
first applied to iteratively estimate the ++ temper- 
ature (T e (0 ++ )) and density (n e ) of the nebula from 
the [O III] and [S II] line ratios, respectively. If the 
measured line ratios correspond to unphysical condi- 
tions (outside the range for which temden is calibrated, 
500 < T e (0 ++ ) < 10 5 K and 1 < n e < 10 8 cm" 3 ), we 
do not calculate the direct abundance. The + temper- 
ature i s then es t imate d using the linear empirical rela- 
tion of lGarnettj (|1992f) . The O II and O III abundances 
are then estimated using the density, ionic temperatures, 
and [O II] and [O III] line ratios follow ing the ionizatio n 
correction factor (ICF) prescription of iShi et al.l (|2006[ ). 
The total oxygen abundance is taken to be the sum of 
these two ionic abundances. 

We employ several strong line diagnostics chosen 
to represent each of the major classes calibrated in 
the literature: R 2 z, N2Q2, N2 , 03N2, and P (see 
lLopez-Sanchez fc Estebanll2010bl . for a recent review). 
First, we apply the R23 oxygen abundance calibration of 
iZaritskv et all (|1994f) . an average of three earlier meth- 
ods, hereafter referred to as "Z94." Z94 is only cal- 
ibrated for the highcr-mctallicity upper branch of the 
i?23-abundance degeneracy (8.4 <log(0/H)+12 < 9.6); 
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TABLE 3 

Nebular Emission Line Fluxes From SN Host Galaxies 



SN 


A v 


[O II] 
A 3727 


[O III] 
A 4363 


H/3 
A 4861 


[O III] 

A 4959 


[O III] 
A 5007 


[N II] 
A 6548 


Ha 
A 6562 


[N II] 
A 6584 


[S II] 
A 6717 


[S II] 
A 6731 


SN Type: lib 
2007br 


170 ± 30 








30 ± 10 




100 ± 30 


38 ±8 


45 ±8 


34 ±7 


2007ca 




350 ± 6 




100 ±2 


84 ± 5 


250 ± 10 


17 ± 1 


360 ± 20 


51 ±3 


54 ±3 


36 ±2 


2010am 


o-o±g : g 


373 ± 3 




100 ±2 


104 ±7 


290 ± 20 


11 ±2 


220 ± 20 


35 ±3 


49 ±4 


32 ±3 


2010cn 


o o +0 - 2 
u - u -o.o 


300 ± 10 




100 ±7 


104 ± 10 


290 ± 30 




270 ± 30 


27 ±5 


45 ±6 


23 ±5 


2011D 


5 6+ 11 

°-°-i.o 






100 ± 30 


60 ±20 


200 ± 50 


120 ± 30 


1500 ± 400 


130 ± 40 


220 ± 70 


140 ± 40 


2011V 


o o +0 - 6 


140 ± 20 




100 ± 40 






80 ±20 


300 ± 100 


100 ± 20 


80 ±20 


80 ±20 


PTF09dxv 5.0±Q7 






100 ± 40 






140 ± 30 


1200 ± 400 


430 ± 60 


320 ± 50 


170 ± 30 



SN Type: 

2002ex 

2004cf 

2005hm 

2006jo 

20061c 

2007az 

2007gg 

2008gc 

2009hu 

2009jf 

LSQ11JW 

PTF09dfk 



lb 

1.3 
0.0 
0.2 
3.1 
0.0 
0.7 
0.4 

5.1 
0.0 
3.5 



+0.5 
0.6 
+0.5 
0.0 
+0.4 
0.2 
+0.3 
0.3 
+0.0 
0.0 
+0.6 
0.6 
+0.3 
0.3 

+0.7 
-0.6 
+0.1 
0.0 
+0.5 
0.5 









5.8 ±0.2 


20.4 ±0.3 


7.2 ±0.6 


100 ± 10 


22 ±2 


25 ±2 


18 ± 1 






100 ± 10 


110 ±10 


300 ± 20 




410 ± 70 


59 ±8 




150 ± 20 




20 ±6 


100 ±4 


180 ± 20 


590 ± 50 




290 ± 50 


27 ±8 


74 ±8 


47 ±7 


450 ± 10 




100 ±5 


110 ±10 


200 ± 20 


55 ± 6 


310 ±40 


190 ± 20 


56 ±6 


46 ±5 






100 ± 20 






60 ± 10 


700 ± 100 


200 ± 20 


76 ±6 


58 ±5 


215 ± 3 




100 ±2 


126 ±7 


400 ± 20 


6±1 


260 ± 20 


12 ±2 


35 ±3 


17 ±2 


310 ± 10 




100 ± 20 


100 ± 20 


240 ± 30 




350 ± 50 


80 ± 10 


120 ± 20 


70 ± 10 


290 ± 20 




100 ±5 


80 ±7 


220 ± 20 


26 ± 3 


320 ± 30 


30 ±3 


52 ±4 


37 ±4 


12 ±5 










22 ± 7 


100 ± 20 


34 ±9 










100 ± 10 


150 ± 20 


190 ± 20 


160 ± 30 


1300 ± 200 


310 ± 40 


110 ±20 


140 ± 30 


300 ± 8 




100 ±8 


50 ±6 


150 ± 10 




260 ± 30 


44 ±6 


40 ±6 




79 ±6 




100 ±4 


110 ±10 


350 ± 30 


30 ±4 


800 ± 100 


110 ± 10 


170 ± 20 


130 ± 10 



SN Type: Ic 








2003cv 










2004ib 


2> j 


630 ± 20 


100 ± 6 


73 ±7 


2006ip 


42 ±2 




2.4 ±0.8 


2006ir 




190 ± 10 


100 ± 3 


79 ±6 


2006tq 


160 ± 10 




50 ± 10 


2007db 


9+ - 4 
u - y -0 .4 


600 ± 20 


100 ± 5 


76 ±7 


2007flF 


30 ±10 






2007gl 










2007hb 


2 1+ - 4 




100 ± 20 


17 ±4 


2007hl 


n c+0.2 
z '°-0.2 


460 ± 10 


100 ± 5 


41 ±4 


2007hn 


140 ± 40 






2008ao 


i =+0.3 
L -°-0.3 


97 ±1 


100 ± 3 


12 ± 1 


2008ik 


1 q+0-6 
■■■■ a -0.5 


370 ± 10 


100 ± 20 




2009nl 


n O + 0.8 

3.8_ .7 




100 ± 20 




2010Q 


0+ 01 
n 7 +0.5 
U - '-0.5 


212 ± 5 


10 ±3 100 ±4 


160 ± 10 


2011it 


260 ± 10 


100 ± 9 


33 ±6 


PTF09iqd 


140 ± 20 






PTF09q 


6 +0 - 3 


10 ±3 






PTFlObip 


393 ± 5 


100 ±4 


85 ±7 


PTFlOvgv 


35.5 ±0.8 




8±1 


PTFllhyg 1.2±°- 6 e 


144 ± 10 


100 ± 30 





30 ±7 


15 ±5 


100 ± 30 


62 ± 7 


36 ±5 


22 ±5 


230 ± 10 


25 ± 5 


570 ± 40 


108 ±8 


128 ±9 


93 ±6 


3.9 ±0.9 


9.8 ±0.9 


100 ± 10 


32 ± 2 


17 ± 1 


14 ± 1 


260 ± 20 


12 ± 3 


360 ± 20 


39 ±3 


54 ±4 


34 ±3 


100 ± 20 




100 ± 20 


80 ±30 




90 ±30 


220 ± 10 


34 ±6 


370 ± 40 


63 ±6 


92 ±7 


66 ±6 






100 ± 20 


63 ±8 


13 ±3 


13 ±3 




20 ± 10 


100 ± 100 


90 ±20 


30 ± 10 


40 ± 10 


17 ±4 


44 ±9 


500 ± 100 


160 ± 20 






91 ± 5 


59 ±4 


600 ± 40 


159 ±9 


137 ±7 


93 ±5 




33 ± 8 


100 ± 40 


91 ±9 


57 ±6 


36 ±5 


37 ± 2 


50 ±4 


450 ± 40 


150 ± 10 


74 ±6 


58 ±5 




90 ± 10 


500 ± 100 


370 ± 30 


170 ± 20 


130 ± 10 






900 ± 200 


120 ± 40 


150 ± 30 




480 ± 40 


10 ±2 


260 ± 30 


12 ±3 


55 ±5 


18 ±2 


110 ± 10 


45 ± 6 


350 ± 50 


100 ± 10 


96 ± 10 


69 ±8 




40 ± 10 


100 ± 30 


50 ±20 


70 ±20 


40 ± 10 




10 ±2 


100 ± 10 


21 ±3 






250 ± 20 


18 ± 3 


340 ± 30 


34 ±4 


72 ±6 


46 ±5 


22 ±2 


9±1 


100 ± 10 


30 ±3 


24 ±2 


19 ±2 


50 ± 10 


80 ± 10 


400 ± 100 


200 ± 20 


70 ± 10 


51 ±9 



SN Type: Ic-BL 

2005nb l.Otoj 

2006nx 0.1±° l 600 ± 20 

20071 • ■ ■ 130 ± 10 

2007ce 0.Qt l 119 ± 1 

2008iu 0.0l %° 110 ± 7 

2010ah 0.7±p;| 380 ± 20 

2010ay O.ol^o 233 ± 7 



100 ± 10 50 ± 10 

100 ±4 120 ±10 
7±3 

9.8 ±0.8 100.0 ± 1.0 200 ±10 

23 ± 4 100 ± 2 210 ± 20 

100 ± 10 60 ± 10 

4.6 ±0.6 100.0 ±0.6 170 ±20 



140 ± 20 • • ■ 380 ± 40 

350 ±40 17 ± 5 290 ± 60 

29 ± 4 • • ■ 100 ± 10 

550 ± 20 2.5 ± 0.4 270 ± 30 

640 ±70 13 ± 2 230 ± 40 

180 ± 20 • • ■ 350 ± 40 

582 ±4 7.7 ±0.6 250 ± 30 



70 ± 10 


70 ± 10 


50 ± 10 


38 ± 7 


80 ± 10 


33 ±6 


23 ±3 


41 ±4 


20 ±3 


7.4 ±0.7 


14 ± 1 


12 ± 1 


11 ± 1 


12 ±2 




33 ±9 


70 ± 10 


43 ± 10 


25 ± 2 


29 ±3 


22 ±2 



SN Type: Undetermined Ibc 

1991R 1-6±8:| 560 ± 30 • ■ ■ 100 ± 6 49 ± 6 

2011gh 2.0t° 4 193 ± 7 ■ ■ ■ 100 ± 10 12 ± 4 

2011hw ■■■ 40 ±10 



180 ±10 41 ± 5 460 ±30 106 ± 8 77 ± 6 54 ± 6 
48 ±6 56 ±7 520 ±80 190 ± 20 93 ± 9 72 ± 7 
23 ± 7 100 ± 20 61 ± 9 40 ± 10 50 ± 10 



Note. — Ay has been derived from the Balmcr decrement as described in the test, but fluxes reported here have not been dcrcddcncd. Fluxes are 
reported relative to H/3 — 100 when possible; Ha or oxygen lines are used for normalization when necessary. 
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TABLE 4 

Metallicities and Ages of SN Host Environments 



SN Host Galaxy 






log(0/H)+12 






Age a 




Direct 


Z94 


KD02 


PP04N2 


PP04O3N2 


PT05 


(Myr) 


CAT TT, irus . Jlh 

oi\ lype. no 
















9007hr 
ZUU I L)I 








8.66 ± 0.09 








9007nn 




8 Af\ 4- OA 
o.*±\j ztz u.u^t 


8.55 ±0.07 


8.42 ±0.02 


o.oo m u.ui 




^94-04 


901 Darn 

£l\J J. W (1111 




Q 47 _i_ ri no 


8.58 ± 0.03 


8.45 ± 0.03 


o.oo zc u.u^ 


SIR -1- 09 
O . J-U zc u.uz 


6.0 ± 0.5 


901 Drn 
iUXUCJl 




e Kyi _(_ n 07 


8.54 ± 0.08 


8.33 ±0.05 


a ok 4. n o^ 
O.ZU I u.uo 


c o« _l n OR 
o. zu ztz u.uu 


fi ^ 4- A 

O.O ^ U.l 


9ni 1 n 








8.30 ±0.10 


O.OO ZH U.UD 




^ A 4- O ^ 

o.^t ztz u.o 


2011V 






9.04 ±0.07 


8.6 ±0.1 






S A 4- Q 
o.4: ^ u.y 


l ± jt uyuxv 








8.64 ±0.09 






q q j-no 

o.o ztz u.y 


SN Type: lb 
















2002ex 








8.52 ±0.04 








zUU4ci 








8.42 ±0.05 


8.31 ± O.Oo 




8.1 ± 0.6 


ZUUOllIll 








8.32 ± 0.09 


o. J-U m u.uo 




7 Q _|_ A 
/ .o ztz U.I 


900fiir\ 
ZUUUJO 






8.93 ±0.05 


8.79 ± 0.04 


o. o t m u.uz 




& A 4- Q 
o.4: ^ u.y 


900filr 








8.59 ±0.05 






8.2 ± 0.8 


9007ar? 

zuu / az 




s ^o 4- o o^ 

O.OU It U.UO 




8.15 ±0.04 


O. ±Z ZJZ u.uz 


oQj.no 

O.O ZtZ U.Z 


^ A 4- O 9 
O.^t ZtZ \J.L 


9OD70-O- 

ZUU i gg 






8.7 ±0.1 


8.52 ±0.05 


s 4n 4- n D 1 ^ 

o.iu zc u.uo 




9.0 ± 0.8 


ZUUOgC 




ft R9 -1- 07 


8.51 ±0.08 


8.31 ±0.03 


£ 9Q -1- fl 09 

o.zy x u.uz 


H 9^ -1- O OS 
O.ZO ZLZ U.UO 


c n -i- n '? 

D.U ZtZ U.O 


900Qhn 

zuuyiiu 








8.64 ± 0.09 








zuuyji 








8.55 ± 0.06 


o.*± / u.uo 




Q 4- S 
y.u ztz u.o 


T ^Ol 1 TW 




8 78 + n Of! 


8.71 ±0.04 


8.46 ± 0.04 


o.*io m u.uo 


c Qfi 4- 0^ 
O.OU ztz u.uu 




r If uyclIK 




8 ^S 4- OS 
o.oo nz u.uo 


8.7 ± 0.1 


8.42 ±0.04 


a 4- n 09 
O.OU zn u.uz 


S "39 4- O OS 

o.oz ztz u.uo 


Q O 4- O R 

y.u ztz u.o 


SN Type: Ic 
















ZUUOfc! V 








8.78 ±0.08 








900Aih 








8.49 ± 0.03 


£ 40 4- 09 
o.iu zc u.u^ 




8.7 ± 0.7 


ZUUOlp 








8.62 ±0.04 








900fiir 




o.uy m u.uj 


8.66 ±0.06 


8.35 ± 0.03 


q on 4. 09 
o . zc u.u^ 


O.O i ^ u.uu 


5.7 ± 0.3 


onnfitn 








8.9 ±0.1 








900 7H h 








8.46 ± 0.04 


Q qq 4. 09 
o.oo zc u.u^ 




7.9 ± 0.6 


90fl7fF 
ZUU i n 








8.78 ±0.07 








2007gl 








8.9 ±0.2 








2007hb 








8.61 ±0.06 


8.82 ± 0.04 




6.3 ±0.7 


2007hl 






8.43 ±0.09 


8.57 ±0.02 


8.57 ±0.01 




8.2 ±0.8 


zUuYrm 








8.88 ±0.10 








OfinSan 
ZUuOaU 




q i q _|_ n 09 
y. j-O m u.uz 


9.02 ± 0.04 


8.63 ± 0.03 


Q _L fl 01 
O. 1 O ZtZ U.Ul 


& AA 4- OA 


fi R 4- 7 
u.u ztz u. / 


ZUUOIK 






8.87 ± 0.09 


8.83 ±0.06 






o.^t zc u.y 


zuuyrn 








8.4 ±0.1 






1 _|_ n c; 

o.± ztz u.o 


901 no 

ZUIUV^ 








8.14 ±0.07 


c HQ -1- OA 

o.uy z\z u.U't 


S 9 4- O 1 
0. z ztz U. 1 


0.0 ztz u.z 


201 lit 




q co _i_ n OQ 
o.oz m u.uy 


8.85 ±0.08 


8.60 ±0.04 


c c;c; _L n no 
O.OO ZtZ u.uz 


c o« _l n 1 
o. zu ztz U. 1U 


1 _i_ n a 
0.1 ztz u.o 


PTFOQinH 








8.7 ±0.1 








PTF09q 








8.52 ± 0.05 








PTFlObip 








8.33 ±0.04 


8.29 ± 0.02 




6.9 ±0.4 


PTFlOvgv 








8.61 ±0.04 








PTFllhyg 






9.03 ±0.08 


8.72 ±0.07 


o.7z ztz U.U4 




7.9 ± 0.9 


SN Type: Ic-BL 
















2005nb 








8.46 ± 0.05 


8.44 ±0.03 




6.4 ±0.5 


2006nx 








8.39 ±0.07 


8.28 ± 0.04 




6.7±0.5 


20071 








8.53 ± 0.05 








2007cc 


7.92 ±0.05 






8.01 ± 0.03 


7.99 ±0.02 


8.3 ±0.2 


3.8 ±0.1 


2008iu 






8.58 ±0.05 


8.14 ±0.05 


8.05 ± 0.03 


8.27 ±0.04 


4.7 ±0.2 


2010ah 




8.6 ± 0.1 




8.31 ± 0.08 


8.32 ± 0.05 




8.0 ±0.4 


2010ay 


8.35 ±0.07 




8.62 ±0.03 


8.33 ±0.04 


8.16 ±0.02 


8.18 ± 0.01 




SN Type: Undetermined Ibc 














1991R 






8.4 ±0.1 


8.54 ±0.02 


8.46 ± 0.02 




8.1 ±0.8 


201 lgh 




8.85 ±0.10 


8.85 ± 0.06 


8.65 ± 0.04 


8.70 ± 0.02 






201 lhw 








8.78 ± 0.07 









Note. — The oxygen abundance, lo g(Q/H) + 12, is estima ted usin g the following diagnosti cs, as described in the text: 
Direct (based on [O i n] A4363), Z 94 JZaritsky et al.|p99l) , KD02 JKewley fe Dopita|[200^ , PP04N2 and PP04O3N2 
iPettini fc Pagellf200l) . and PT05 dPilvugin Thuanll200g) . 

a The age of the young stellar population has been estimated based on the rest-frame equivalent width of the Hf3 line, 
following Lcvcsquc et alj H2010al) ( sec Scction|3.5|l. 
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in every case where we apply Z94, the [N II] /[O II] 
ratio suggests an upper branch solution. Second, 
we apply t he [N II] /[O II] oxyge n abundance cal- 
ibration of IKewlev fc Dopital (|2002t ) (as updated by 
Kewlev fc Ellison|2008D hereafter referred to as "KD02." 
Kewlev &: Dopital (|2002D synthesize a variety of modern 
photoionization models and observational calibrations to 
produce recommendations for producing an abundance 
estimate given different permutations of available emis- 
sion lines and uses the [N II] /[O II] ratio to break the 
degeneracy between the upper and lower branches of 
R 23 . Third, we apply the empirical [N II] /Ha ("N2") 
and [O III] /[N II] ("Q3N2 " ) oxy gen abundance cal- 
ibrations of iPettini fc Pagell (120041 ). hereafter referred 
to as "PP04." Fourth, we apply the excitation pa- 
ra meter ( "P method " ) oxy gen abundance calibration 
of iPilvugin fc Thuanl (|2005| ). hereafter referred to as 
"PT05." P is c alculated from the ratio of [O III] to 
([O II] +[0 III]) (|Pilvuginll200l . and the [N II] /[O II] 
ratio is used to break the R23 degeneracy. 

There are well-known offsets between the diagnostics 
which are particularly large between empiri cally and 
theor etically-calibrated diagnostics (see e.g. iStasinska 
2002). However, the relative metallicity difference mea- 
sured between a given pair of galaxies in different di- 
agnostics is consistent with an rms scatter typically ~ 
0.07 de x, and 0.15 dex between t he most discrepant diag- 
nostics (jKewlev fc Ellisonll2008f ). Comparing the metal- 
licities we measure for the same host galaxies in different 
diagnostics, we find discrepancies consi s tent w ith the rms 
scatter reported by IKewlev k, Ellison] (|2008l ). Hereafter 
we refer to this uncertainty intrinsic to the diagnostics 
as the "systematic uncertainty;" we do not factor the 
systematic uncertainty into the metallicities reported in 
Table IH but we do consider this systematic uncertainty 
in our statistical analysis (Section I3.4|) . The system- 
atic uncertainty is typically as large as the "statistical 
uncertainty" associated with the line flux and Ay mea- 
surement errors; for example, the median statistical un- 
certainty associated with our PP04N2 measurements is 
0.09 dex. 

The galaxy sample varies with the diagnostic cho- 
sen. The PP04N2 diagnostic can be applied to nearly 
all galaxies in our sample (N = 50), while the sim- 
ilar PP04O3N2 diagnostic can only be applied to the 
31 galaxies which have measurements of the potentially 
fainter [O III] and H/3 lines. The [PT05,Z94] diagnostics 
can only be applied to [13,12] galaxies, due to their strict 
dependence on the full complement of [O II] and [O III] 
lines. We can apply the direct diagnostic to only 2 galax- 
ies due to its reliance on the auroral line which is weak 
at the metallicity regime probed here. However, the N2 
diagnostic produces relative metallicit y estimates consis- 
tent w ith the other strong line methods IKewlev fc Ellisonl 
(|2008t ). N2 is also less sensitive to systematic effects: it 
employs lines with a very small wavelength separation 
(so extinction correction may be neglected) and the ef- 
fect of underlying absorption is less important (because 
the absorption equivalent width of Ha is typically equal 
to o r less than H/3, despite t he ~ 3 times stronger flux, 
e.g. iBrinchrnann et al"1l2004D . However, in N2 diagnos- 
tics no correction is made for ionization parameter and 
the N2 ratio can saturate at high met allicities when ni- 
trogen becomes the dominant coolant (|Kewlev fc Dopital 



12001 . 

3.3. Metallicity Distribution of SN Ibc Progenitor 
Environments 

Using the PP04N2 abundance diagnostic, we measure 
metallicities for 12 SN lb host galaxies and find values 
ranging from log(0/H)+12 = 8.15 — 8.79, with a median 
abundance and standard deviation of log(0/H)+12 = 
8.48 and 0.16 dex, respectively (Figure 3]). The char- 
acteristics of the 21 SN Ic host galaxies in this di- 
agnostic are similar, with metallicities ranging from 
log(0/H)+12 = 8.14 - 8.88 and median and standard 
deviation of log(0/H)+12 = 8.61 and 0.22 dex. Among 
the 7 SNc lib for which wc measure PP04N2 metallici- 
ties, we find a range of log(0/H)+12 — 8.30 — 8.66 with 
median log(0/H)+12 = 8.46 and the standard devia- 
tion is 0.18 dex dex, similar to the SN lb. In contrast, 
the 7 SN Ic-BL host environments typically have lower 
metallicities, with median log(0/H)+12 = 8.34 and stan- 
dard deviation of 0.21 dex. The minimum metallicity for 
SNe Ic-BL host galaxies (log(0/H)+12 = 8.01) is only 
^0.1 dex lower that of SN lb and Ic host galaxies, but 
the highest metallicity measured for an SN Ic-BL host 
galaxy (log(0/H)+12 = 8.53) is similar to the median 
for SN lb and Ic. 

We report the results f rom the other di agnostics in 
Table El Wc also note that lStoll et~aTI (|2012D have devel- 
oped an emperical, linear conversion between oxygen and 
iron abundance. Applying this to our host environment 
metallicities yields [Fe/H] = [-0.6±0.2, -0.4±0.3, -0.6± 
0.2, -0.8 ±0.3] using our full sample of SNe [Ib,Ic,IIb,Ic- 
BL]. 

Table[5]illustratcs that the effect of restricting our sam- 
ple to only those objects with Gold classifications and ex- 
plosion site spectroscopy is modest. Among the 21 SNe 
Ic host galaxies in our sample for which the PP04N2 di- 
agnostic was applied, 13 have Gold classifications and 
explosion site spectroscopy. The difference between the 
median metallicity measured between this subsamplc and 
the full sample is -0.01 dex. For SNe [Ib,IIb,Ic-BL], the 
subsamplc fraction is [8/12, 5/7, 3/7] and the difference in 
the median is similarly small at [—0.05, —0.04, 0.02] dex. 
These differences, smaller than the systematic uncer- 
tainty associated with the abundance diagnostics (Sec- 
tion 13. 2p , generally reflect the lower metallicities of the 
SN explosion sites as compared to the galaxy nuclei. 

Finally, there are several SN host environments at 
which we performed spectroscopy, but could not measure 
metallicities. For 7 objects, we detect narrow emission 
lines in the host galaxy, but they are not sufficient to 
estimate the metallicity in any of the strong line diag- 
nostics: 2002gz (lib), 2008fi (lib), 2008im (lb), 201 lbv 
(lib), 2011cs (Ic), 2011ip (Ic), and PTFlOaavz (Ic-BL). 
For 2 objects, we did not detect any narrow emission lines 
in our host galaxy spectrum: 2004ai (lb) and 20101z (Ic). 
With the exception of SN 2004ai, these host galaxies are 
not at exceptionally high redshifts with respect to the re- 
mainder of the sample, and poor S/N in the spectroscopy 
is due to the intrinsic low luminosity of the host galaxies. 

We place upper limits on the host environment metal- 
licity of these SNe using archival photometry co mbined 
with the L - Z relation of iTremonti et all (|2004r f rr l and, 

11 Converting the Trcmonti et al. (2004) L — Z relation from 




Fig. 4. — Cumulative metallicity functions for SN host galaxies from this study, using the PP04N2 diagnostic. The 16th and 84th 
percentiles (lcr limits) of the probability distribution for the CDFs are illustrated by the dashed lines, as computed by propagation of 
the uncertainties in the line flux measurements and the 0.07 dex diagnostic systematic error. For compa rison, we show th e distribution 
of PP04N2 explosion site mctallici ties of SN II discovered b y the untargeted PTF surv ey as measured bylStoll et"a l. (2012), and that of 
low-redshift (z < 0.3) LGRBs from lLevesque et al.l H2010al) , iLevesque et all 12010d 1. and IChornock et alj 1 120101) . 



when Ha is detected, the the 3<r upper limit measured 
for the [N II] flux combined with the PP04N2 diagnostic. 
We summarize this investigation in Table [5] and note a 
few special cases here, but we do not use these metal- 
licity limits in our figures or statistics except where ex- 
plicitly noted. SN 2008im occurred ~ 8 kpc from the 
nucleus of the Sb galaxy UGC 02906 (jOksanenl [20081 
with z = 0.008). In our spectrum, light from an older 
stellar population dominates over any signature of star 
formation. Given that the host galaxy has an absolute 
magnitude Mb ~ — 19 mag, the explosion site is likely 
to be of sub-solar metallicity. The only one of these nine 
host environments for which we cannot place useful con- 
straints on metallicity is that of SN 2004ai, reported as a 

the T04 to the PP 04N2 scale and adopting log(Q/ H) w + 12 = 8.9 
on the T04 scale (Dclahayc & Pinsonncault 2006), a galaxy with 
Mb ~ —19.8 mag should have solar metallicity. 



SN Ic bv lRiello et all (I2004al) . but revised to SN lb (M. 
T. Botticella, private communication). The supernova 
spectrum indicates z ~ 0.59, making it by far the most 
distant SN in our sample. Ha is redshifted out of our 
spectral range and nothing is visible in DSS images at 
this position, from which we infer that Mb > —23 mag, 
which does not allow us to distinguish between sub- or 
super-solar mctallicities. The remaining seven host envi- 
ronments are constrained to be at sub-solar mctallicities, 
although for SN 2008im and 20101z we do not have suf- 
ficient S/N in Ha to place spectroscopic limits on the 
metallicity. 

3.4. Statistical tests on metallicity distributions 

We apply the KS test to our metallicity measurements 
(Table [7| and interpret pks < 0.05 to indicate statis- 
tically significant evidence for a difference in the par- 
ent populations of the two sets being compared. When 
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TABLE 5 

Statistics of Measured Host Galaxy Metallicity (log(0/H)+12) Distribution 



Diagnostic 



Minimum 



Maximum 



Median 



SNe lib 
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Note. — The statistical properties of the distribution of mctallicitics measured for the host 
galaxies, divided by SN type. The first values listed represent all the SNe in the sample; the 
values in parenthesis reflect only those SNe with secure typing (Gold sample; Section 12. Il l 
and explosion site spectroscopy (Section \2.2t . The oxygen abundance diagnostics applied 
are described in Section 13.21 SNe of type "Ib/c" have uncertain typing. The medians and 
standard deviations (cr) have been calculated from the MCMC samplings. 



0.20 

0.08 



Monte Carlo simulations^ indicate that pks rises above 
this threshold (> 0.05) in at least 14% of trials, we refer 
to this as "marginal" evidence of statistical significance. 

We find no significant difference between the metallic 
ity distribution of SN lb and Ic, with j»ks = O.lOjl 
using the PP04N2 diagnostic (p K s = 0.07^^5 if re- 
stricted to Gold classifications and explosion site metal- 
licity measurements ) . Th is contrasts with the find- 
ing of IModiaz et al.l (|2011[ ) , that the distributions dis- 
agree at the pks — 0.01 level using the equivalent 
PP04 O3N2 diagnostic, but is c onsistent with the null re - 
sult of I Anderson et al.1 (f2010h : iKeTTv fc Kirshnerl (I2TIT1 : 
iLeloudas et al.l (1201 1[ ). We note that whe n we apply the 
KS test to the sample of SN lb and Ic from lModiaz et aLl 
(2011) using our methodology, we find pks = 0.08lQQg. 



The difference is that we separate SN lib and Ibn from 
lb, which changes the metallicity distribution and re- 
duces the SN lb sample size by 5 objects. 

We find different metallicity distributions for SN Ic and 
Ic-BL, with pks = O.OIIq 1 ° l using the PP04N2 diagnos- 
tic. Given the smaller sample size and lower median 
metallicity of the SN lb, the evidence for a difference in 
the Ic-BL and lb metallicity distributions is not signifi- 
cant (pks = 0.21^Q ^g). We find no evidence for a signif- 
icant difference between the SN lb and lib populations 



KS 



0.75 



+0.22 

-0.23 



)■ 



3.5. 



12 Following ILeloudas et aTl l|2011l ). we incorporate the uncer- 
tainty in the individual metallicity estimates (but not the diag- 
nostic systematic uncertainty) by repeating the KS test through 
Monte Carlo simulations where we sample from the full probabil- 
ity distribution for the metallicity of each host galaxy. 



Young stellar population ages 

We estimate the a ge of the young stellar population 
using the method of iLevesque et all ()2010al ). assuming 
an insta ntaneous-burst star formatio n history (for a re- 
view see iStasihska fc Leitherer| [l996') . The age estimate 
is based on the rest frame equivalent width of the H/3 
line (Wh/3) and evolutionary synthesis models for star- 
burst galaxies based on the Geneva HIGH evolution- 
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TABLE 6 

Host environments without metallicity measurements 



SN 


SN Type 


z 


M B (mag) 


Z 


(phot.) a 


Z 


(spec.) b 


2002gz 


lib-pec (G) 


0.085 


> 


-18 


< 


8.4 


5 8.5 


2004ai 


lb (S) 


0.590 


> 


-23 


< 


8.9 






2008fi 


lib (G) 


0.026 


> 


-15 


< 


8.2 


S > 


1.2 


2008im 


lb (G) 


0.008 




-19 




8.5 






20101z 


Ic(G) 


0.090 c 


> 


-12 


< 


8.0 d 






201 lbv 


lib (G) 


0.072 


> 


-17 


< 


8.3 


< 1 


3.3 


2011cs 


Ic (G) 


0.101 


> 


-18 


< 


8.1 


< ! 


i.l 


2011ip 


Ic (G) 


0.051 


> 


-17 


< 


8.3 


S '' 


S.5 


PTFlOaavz 


Ic-BL (G) 


0.063 




16.5 




8.2 


S ' 


3.1 



Note. — The security of the spectral classification (Silver, S, or Gold, G; 
see Section |2.H is indicated in parenthesis. Unless otherwise noted, we retrieve 
photometry for these host galaxies from NED, the NASA/IPAC Extragalactic 
Database operated by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space Administration. 
When photometry is not available in NED (not detected in DSS), we assume 
rriB > 20 mag. Rcdshifts are taken from NED or our own spectroscopy, except 
where noted. 

a Metallicity (log(0/H) +12) limit implied by the L — Z relation using the pho- 
tometry listed, on the PP04 scale. 

a Metallicity (log(0/H) +12) limit implied by the PP04N2 diagnostic using the 
3er upper limit measured for the [N n] flux, stated when Ha is detected. 
c I Drake et"aH J2011a|) 

The absolute magnitude limit places this host galaxy below the range over 
which the L — Z relation of [Trcmonti el al. (2004) is calibrated (~ 8.0 — 9.2 on 
the T04 scale). 



TABLE 7 





KS TEST 


p-VALUES FOR SNE IBC 






lib 


lb 


Ic 


Ic-BL 


lib 

lb 

Ic 

Ic-BL 


67+ ' 27 
u -°'-0.30 

u - iD -0.09 


75+ - 22 
u ' ,o -0.23 

07+ ' 21 
u - u '-0.05 


34+ ' 17 

o in+ - 20 

u ' iu -0.08 


42+ ' 00 
21 +0 ' 31 

o m+ 003 

u ' ui -0.01 



Note. — The values in the table are KS test p-values 
(Pks) f° r the probability that the measured host galaxy 
mctallicities of the SNc of the two indicated types were 
drawn from the same parent populations. Above the diag- 
onal, all objects in our sample arc considered; below the 
diagonal, only objects with Gold spectroscopic classifica- 
tion and explosion site spectroscopy arc considered. The 
PP04N2 diagnostic is used for metallicity estimation. The 
KS tost is performed only if metallicity measurements arc 
available for > 4 SN of each type. The upper and lower 
limits listed are the 16th and 84th percentile values of the 
results of the KS test Monte Carlo simulations. 



ary tracks (jMevnet et al.lll994t ISchaerer fe Vaccalll998f ). 
We use the PP04N2 diagnostic to break the metallicity 
degeneracy, restrict the sample to spectra obtained at 
the explosion site, and exclude objects with supernova 
flux contamination (SNe 2010ay, 2011bv, 2011gh, 2011ip, 
LSQ11JW, PTFlOaav). 

As illustrated in Figure \E\ wc find that the ages of 
young stellar populations in Type lb and Ic SN host 
environments are not discrepant, with median ages of 
[8.1,7.4] Myr and standard deviations of [1.5,1.1] Myr 
for N = [6, 10] SN [Ib,Ic] with explosion site spectra, 
and pks = 0.21. However, the ages of SN Ic-BL environ- 
ments (N = 4) are somewhat lower, with a median of 
5.6 Myr and standard deviation of 1.6 Myr. Given the 
small sample size, the KS test cannot confirm statistical 
significance, pks = 0.34. Similarly, the age of SN lib 
environments seems to be lower than that of SN lb and 
Ic, with median age 6.0 Myr and standard deviation of 
1.2 Myr, but the sample size is small (N = 5). 



We note that the distribution of ages for SN Ic-BL host 
environments in our sample are similar to that of the low- 
redshift (z < 0.3) LGRB hosts studied in lLevesque et all 
(|2010al) . which have a median age of ~ 5 Myr. The 
KS test does not suggest that the LGRB and SN Ic-BL 
age distributions are significantly different (pks = 0.48), 
but with only 5 low- 2; LGRBs the sample sizes are small. 
Moreover, the observation that SN Ic-BL and lib occur in 
y ounger stellar popu l ations is consistent with the finding 
of lKeilv fe Kirshnerl (|2011h that these types of SNe have 
bluer explosion site u' — z' colors. 

With some caveats, we can interpret the young stellar 
population ages in terms of the lifetime o f massive stars. 
For con sistency with the W^p models of iLevesque et al.1 
(2010a), we use the stellar lifetimes from the "Geneva 
HIGH" (high mass loss) solar-m etallicity (Z = 0-02) 
evolutionary tracks computed by iMevnet et al.l (|1994f ) 
(shown in Figure [5]). The measured ages therefore imply 
progenitor stars with initial masses Mj ~ 20 — 40 M Q 
for SN lb and Ic. This is similar to the mass range 
for SN Ib/c progenitors expected from stellar evolution- 
ary theory and indicated by the relat ive rate of SN Ibc 
and SN II £B oissicr fc Prantzod 120091 ). but may be dis- 
crepant with progenitor non-detection up per limits from 
pre-explosion imaging (for a review, see iSmarttl l2009f l . 
For SN lib and Ic-BL, the young stellar population ages 
imply progenitor stars that are somewhat more massive, 
Mi < 60 Mq. While stars of a given mass are longer- 
lived at lower metallicity in the Geneva models, and the 
SN Ic-BL progenitor stars are found in lower metallicity 
environment, the effect is small because the lifetimes vary 
by ~ 10% over a factor of ~ 3 in metallicity. However, 
because the ages are estimated assuming an instanta- 
neous burst of star formation and neglecting any ongoing 
star formation, the possibility of younger progenitor stars 
within the population is not precluded. Further assump- 
tions inherent to this Wh/3 diagnostic include complete 
absorption of ionizing photons, spatially uniform dust 
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SN Type A+10 K+ll L+ll M+H S+12 TW 



Mi (M s ) 
40 25 




lib 

5 

lb 

6 

Ic 

10 

Ic-BL 
4 

LGRB 
5 



TABLE 8 
Host Galaxy Samples by type 



Targeted SNe 

lib 
lb 

Ib/c 
Ic 

Ic-BL 



1 
10 

3 
14 





13 
10 

2 
23 

5 



3 
7 
1 
10 
6 



0.005 0.011 0.016 0.012 



3456789 10 
W HS Age (Myr) 



Fig. 5. — Cumulative young stellar population age distribution 
for the SN Ibc host galaxies, based on the rest frame equivalent 
width of the H/3 line (Wnp) measured from explosion-site spec- 
troscopy only. The color coding and number of objects of each 
type are listed on the right. The low-redshift (z < .3) LGRB host 
environment ages come from Lcvesquc ct al. (2010a). The top axis 
illustrates the initial masses of rotating Wolf-Rayet stars with lifc- 
times equivalent to the stellar population ages, as predicted by 
IJVlevnet etaTI JT994T) . 

extinction in both nebular and stellar emission regions, 
and it is based on evolutionary tracks that do not include 
the effects of stellar rotation. 

3.6. Wolf-Rayet star populations 

Broad Wolf-Rayet (WR) features ( "bumps" ) in galaxy 
spectra, reflecting the existence of evolved, massive stars 
(M > 25 M Q ), can be used to character ize the na- 
ture o f ongoing star formation i n the galaxy (I Allen et al] 
19761: iKunth fc Sargent] fTflSTt ISchaerer fc Vaccal 119981 : 
Schaerer et al.l I1999T ) . In particular, the "blue bump" 
which is primarily due to the He II A4686 line is an in- 
dicator of late-type WN Wolf-Rayet (WNL) stars. A 
visual inspection of the spectra of our SN Ibc host envi- 
ronments does not reveal any recognizable blue bumps. 
We estimate the 3a upper limit of the flux in the blue 
bump feature as 3 times the rms of the continuum flux 
in a 40 A window at its locat ion. 

We f ollow the method of ILopez-Sanchez fc Estebanl 
(|2010a[ ) to place limits on the fraction of WNL stars in 
the young stellar population, WNL/ (WNL + O), based 
on the ratio of the upper limit flux of the blue bump 
to H(3. For WR galaxies in their dataset, this fraction 
ranges from WNL/(WNL + O) - 0.03 - 0.3, smaller 
than we can constrain with most of our spectra due 
to the continuum S/N. However, there are 3 host en- 
vironments for which we can rule out WNL popula- 
tions at that level (2007az, 2008iu, 2010Q; all explo- 
sion site spectra) and one for which we can rule out 
WNL/(WNL + 0) > 0.1 (2007ce). Significantly higher 
S /N spectroscopy or narrow-band imaging could provide 
stricter constraints for typical SN host galaxies. This 
analysis demonstrates that investigations of the WR pop- 
ulations of SN host environments requires a significantly 
different observing strategy than a study designed for 
strong line metallicity measurements, and will likely be 
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15 
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z 


0.005 


0.015 


0.022 


0.017 


0.036 
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Note. — The number of host galaxies in each of the litera- 
ture samples divided by SN type, and median redshifts, z. The 
numbers are given for SNe discovered by untargeted (targeted) 
surveys. SNe of type "Ib/c" have uncerta in typing. The sample s 
come from the following references: A+10, [Anderson" et al.l (2010); 
K+ll, [Kelly <k. Kirshncr 1201 1]) (only ob jects wit h metallicities 
measu red); L +ll, Lcloudas ct al. 1201111 ; M+ll. [Modiaz et all 
120111) ; S+12, ISto"ll et al.N20121) ; TW. this work. 



TABLE 9 

Overlap between SN Ibc host galaxy samples 



A+10 K+ll L+ll M+ll S+12 TW 



A+10 28 7 

K+ll 67 3 12 3 4 

L+ll 20 5 5 

M+ll 47 5 

S+12 9 4 

TW 58 



Note. — The numbers on the diagonal represent the to- 



tal number of SNe Ibc in each sample; the numbers above 
the diagonal represent the intersection of the samples. The 
samples arc the same as in Tabled] 

limited to nearby SN host galaxies. 

4. COMBINED SN Ibc DATASET 

Next we combine our dataset with those of pre- 
yious /concurrent spectr osco pic studies of S Ne Ibc 
(I Anderson et all \2(M. Kelly fc Kirshnerl [2?)TTl 
iLeloudas et al] I2011L iModiaz et al] 120111 and this 
work). We summarize the characteristics of these 
surveys in Table [8] 

4.1. Criteria of the Combined SN Ibc dataset 

For the purposes of assembling a statistical sample, 
we consider only measurements made on the PP04 scale 
(see Section I3.2|P^I . Among these five samples, there 

13 For our observat ions and those of Anderson et al. (20l7j) and 
ILeloud as et al. (2011). we employed measurements made using the 
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are metallicity measurements for 171 unique SN Ibc host 
galaxies, of which 133 had metallicity measurements pre- 
dating our study. To ensure consistency between au- 
thors when propagating uncertainties, we add in quadra- 
ture a representative systematic uncertainty of 0.07 dcx 
(see S ection | 3.2I) to the me tallicit y estimates from ou r 
stu dy. ILeloudas et al l (120111) . an d IModiaz et al.l (120111 ); 
for iKellv fc Kirshnerl (pOlTfT and iStoll et al.l (|2012l ). who 
do not report metallicity uncertainties, we assume this is 
the sole uncertainty. 

Following the auth ors' own evaluations , we consider 
all observat ions from IModiaz et al.l (|2011l ) and all but 
three from ILeloudas et al.l ( 20111) to be at the explo- 
sion site. If we were to apply o ur own criteria (Sec - 
tion !2.2[) . 4 of the observations from lModiaz et al.l (|2011l ) 
would not qualify as explosion site due to their high red- 
shift (SNe 2007 .j y, 200 7qw, 2005kr, and 2006nx). For 
lAnderson et al.l (|2010l ). following Section 12.21 we con- 
sider the 25/28 observations where the spectrum was 
extracted < 2 kp c from the SN pos it ion to be at the ex- 
plosion site. For IKellv fc Kirshnerl (|201lD . who employ 
3" fiber spectroscopy from the SDSS, we consider all of 
the observations to be nuclear rather than explosion site 
measurements, as they sample the global properties of 
the host galaxy. We have investigated the discrepancies 
between those host galaxies whose metallicity measure- 
ments were performed by multiple authors and they are 
typically < 0.15 dex. When SNe have metallicities re- 
ported by multiple authors ( see Table [9]), we adopt the 
average and use only explosion site metallicities where 
possible. We neglect any systematic uncertainty intro- 
duced by differences in instrumental characteristics and 
spectroscopic analysis between the samples. 

Of this combined sample, 75 of the SNe were discovered 
by untargeted searches and 114 of the metallicity mea- 
surements come from explosion site spectroscopy. Before 
our observations, these sample sizes were only 37 and 84, 
respectively. Our results therefore approximately double 
the number of untargctcd SN Ibc for which host galaxy 
metallicity measurements have been published, and sig- 
nificantly increase the number of explosion-site metallic- 
ity measurements. 

Figure |5] compares the metallicity measurements made 
in this work to 18 previous measurements of the same 
host galaxies from the combined sample. Generally, dis- 
crepancies are small, with a mean residual of 0.06 dex 
and rms of 0.15 dex. In the median, the discrep- 
ancy is 1.6 x the statistical error in our measurement. 
The largest outliers can be explained by clear differ- 
ences in spectroscopic methodology between our sample 
and the previous works. For example, for PTF09q we 
measured log(0/H ) + 12 = 8 . 52 in our explosion site 
spectrum, whereas IStoll et aTl (|2012l ) measure a higher 
log(0/H) + 12 = 8.81 from their SDSS spectrum due to 
the nuclear placement of the Sloan fiber (~ 5 kpc from 
the SN site). There are 5 cases where we can compare an 
explosion site measurement from our sample to an equiv- 
alent measurement from the literature, and in those cases 
the rms discrepancy is only half as large (0.08 dex). 




8.2 8.3 8.4 8.5 8.6 8.7 8.8 
log(0/H)+12 (PP04N2, Previous) 



Fig. 6. — Comparison of SN host galaxy metallicity measure- 
ments this work and previous studies in the combined sample. The 
overlap between these samples is described in Table [9] The error 
bars correspond to the statistical error estimated in this work, from 
propagation of the line flux errors. The marker shapes defined in 
the legend as "X, Y" correspond to the spectroscopy type (explo- 
sion site, "E" and nuclear, "N") of the previous work (X) and this 
work (Y). The colors correspond to SN types as in Figure [4] 





lb, N=47 (17) 




lib, N=25 (7) 




Ic, N=68 (21) 




Ic-BL, N=18 (11) 



JO 

CO 
O 




PP04N2 diagnostic; for I Kelly &T Kirsmicr ( 20TE|) and lModiaz et al.l 



1120111 ). who do not report PP04N2 measurements, we instead em- 
ploy the PP 04O3N2 measurements a nd apply the small transfor- 
mation from lKewlev fc Ellison] pOOST ). 



8.0 8.2 8.4 8.6 8.8 9.0 

lopj(0/H)+12 (PP04N2) 

Fig. 7. — Continuous probability distribution of SN Ibc host en- 
vironment metallicities from the combined sample. The curves 
are sums of Gaussians defined by the metallicity measurement and 
their associated uncertainties. Top: All SN Ibc (sample size N 
noted in legend); Bottom: Only SN Ibc from untargeted surveys 
with explosion site metallicity measurements (sample size in paren- 
thesis in legend). 

4.2. SN lb vs. Ic metallicities 

Figure [7] illustrates the metallicity distribution for all 
the SN Ibc host galaxies from the combined sample. First 
we consider the distributions with no selection criteria 
for SN discovery or spectroscopy methodology. With 
sample sizes of [47,68], we find the median metallicity 
of SN [Ib,Ic] to be log(0/H)+12 = [8.56,8.65]. The dif- 
ference between the distributions is not statistically sig- 
nificant (pks = 0.08^0 ofi) an< l tnc difference in the me- 
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Z /3 




.1 8.3 8.5 8.7 8.9 8.1 8.3 8.5 8.7 8.9 

log(0/H)+12 (PP04) 

Fig. 8. — Cumulative mctallicity functions for SN host environments from the combin ed s ample, using the PP04 diagnostic. The colors 
corresponding to SN types are given in the first panel; the sample sizes are given in Table [Tol As labeled, the panels correspond to different 
cuts on the discovery mode of the supernova (galaxy-targeted or untargcted) and the method of spectroscopy (nuclear or explosion site) . 
The dashed colored lines arc la limits for the CDFs as in Figurc[4] but incorporating the systematic uncertainty of 0.07 dex (Scction l3.2H . 
The dashed vertical lines illustrate Zq and Zq/3. 



dians is small relative to the width of the distributions 
([0.16,0.13] dex standard deviation). 

If instead we only consider SNe discovered by tar- 
geted SN searches, we find SN host environments with 
systematically higher metallicities (as quantified in Sec- 
tion 15.11) and we find that the differences between the 
mctallicity distributions for SNe of different types are 
reduced (Figure [8]). The median mctallicity of the 
N = [26, 39] SN [Ib,Ic] from the targeted searches is 

log(0/H)+12 = [8.66,8.65] (p K s = 0.65^). Among 
the subsample with spectroscopy at the SN explosion site 
(N = [19,21] SN [Ib,Ic]) we find median metallicities of 
log(0/H) + 12= [8.62,8.62] ( Pks = 0.63+°;^). 

Looking exclusively to untargeted SNe, for N = [21, 29] 
SN [Ib,Ic] we find median metallicities of log(0/H)+12 = 
[8.46, 8.61] (Figure[8]). While this difference is marginally 
statistically significant (pks = 2.0j^' 7 x 10~ 2 ), it is bi- 



ased by an unequal numbers of galaxy- nucleus versus ex- 
plosion site spectroscopy in the SN lb and Ic samples 
(19 and 28% nuclear spectroscopy, respectively). This 
sample construction bias raises the metallicities of SN Ic 
relative to SN lb (see Section [53]) . If the sample is re- 
stricted to only explosion site measurements, the median 
difference is similar (~ 0.15 dex), but the difference in the 

full distribution is not significant (pks = 0-06i 'u5)- The 
difference in the distributions is most apparent at the 
high-mctallicity end, where very few SN lb discovered 
by untargeted surveys are found at super-solar metallici- 
ties. However, the explosion site spectroscopy in Figured] 
illustrates that SN lb discovered by targeted searches do 
occur in super-solar metallicity environments, and there- 
fore their absence among the untargeted objects must be 
attributed to small sample size. 
How many observations would be required to distin- 
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Fig. 9. — Simulations demonstrating the sample size N of SN Ibc 
necessary to distinguish a difference in the median metallicity of 
SN lb and Ic of [0.1,0.2] dex ([red, blue]) assuming an underlying 
Gaussian distribution with a = 0.2 dex. The shaded regions illus- 
trate the la range in j>ks over random draws from the distribution. 
The horizontal dashed line illustrates the threshold for statistical 
significant, pks = 0.05 at the la level. The vertical lines mark the 
sample size necessary to distinguish the difference in the metallicity 
distribution to this threshold. 

guish between a true difference in the underlying distri- 
bution of progenitor metallicities for SN lb and Ic? If we 
assume that the metallicity distribution of SN Ic progen- 
itors is enriched by 0.1 dex with respect to SN lb progen- 
itors (as we find in Section [373]) and that both distribu- 
tions are Gaussians with standard deviation a = 0.2 dex, 
then we can randomly sample from these distributions to 
investigate the value of pas we would infer from studies 
of different sizes (sec Figure [5]). We take statistical sig- 
nificance to be indicated by the KS test when pxs < 0.05 
and we assume a sample ratio Ni c /Nib = 1.6. 

We find that, in the absence of systematics, only 
N ~ 50 (20) SN Ibc would be required to distinguish 
the discrepancy in 85% (50%) of trials given a dif- 
ference i n the median meta llicity of 0.2 dex (as sug- 
gested bv lModiaz et al.ll201lD . The combined sample in- 
cludes N = [47, 68] SN [Ib,Ic] metallicity measurements; 
N = [21, 29] of these come from untargeted SN searches 
and therefore have substantially reduced systematics (see 
Section [OJ- Therefore, the observations to date should 
be sufficient to distinguish a median metallicity difference 
of 0.2 dex between SN lb and Ic (even among exclusivcly- 
untargcted surveys), which is not supported by the data 
in the combined sample. However, a much larger sample 
of TV ~ 200 (100) observations of SN lb and Ic would be 
necessary to distinguish a median difference of 0.1 dex in 
the progenitor distribution in 84% (50%) of trials. There- 
fore a sample > 2x as large as the combined sample 
would be required to unambiguously distinguish a rela- 
tively subtle discrepancy of 0.1 dex. 

The host environment metallicity measurement s made 
for 34 untargeted Type II SNe by iStoll et all (pOll ) 
(shown in Figure U) constitute an interesting proxy for 
the metallicity distribution of the massive star progen- 
itors of core-collapse SNe. Applying the KS test using 
all the untargeted objects in the combined sample, we 



do not find a significant difference between the metallic- 
ities of SN Ic and SN II (p KS = 0.52), although there 
is a significant difference between SN II and [Ib,Ic-BL] 
(pks = [1 x 10~ 3 ,8 x 10~ 4 ]). Combining all the SN Ibc 
in our sample, we find that they have a median metal- 
licity 0.13 dex lower than the SN II and the difference 
in the full distribution is significant at pks = 0.01. A 
higher median metallicity for SN II than Ib c would con- 
flict with several previous findings (e.g. | Prieto et al.l 
20081 iBoissier fc Prantzoi [200l lAnderson et al.l 120101: 
Kelly fc Kirshnerl 1201 It ) and would not be consistent 
with metal-line driven winds stripping the progenitors 
of SN Ibc, unless combined with a significantly more 
bottom-heavy IMF at low metallicities. This discrepancy 
may instead indicate differences in sample construction. 
The relatively bright SN Ic-BL may be over-represented 
in the combined sample relative to a volume limited 
survey (see Section 15. 4[) , biasing the SN Ibc metallicity 
downward, and some low-luminosity SN II and Ibc host 
galaxies do not have metallicity meas urements (see Se c- 
tion 15.51 in this work and Section 3 of IStoll et al.H2012T) . 

4.3. SN Ic vs. Ic-BL metallicities 

When no cut is placed on SN discovery or spectroscopy 
characteristics, the SN Ic-BL in the combined sample 
(N = 18) have an appreciably smaller (> 0.2 dex) 
median metallicity (log(0/H)+12 = 8.36) than other 
SN Ibc, and the difference relative to the SN Ic distri- 
bution is significant with j?ks = 1 -3 jlf x 10 -4 (see Fig- 
ure [H]). Restricting our scope to SNe discovered by tar- 
geted searches, the sample of targeted SN Ic-BL is limited 
(N = 6). The SN Ic-BL again show a median metallic- 
ity (log(0/H)+12 = [8.45]) ~ 0.2 dex smaller than other 
SN Ibc, but the difference compared to the SN Ic has only 
marginal statistical significance (pks = 4.4^ 8 x 10~ 2 ). 
However, for A = 12 SN Ic-B L from untargeted su rveys 
(primarily from this work and lModiaz et all (|2011l) ). we 
find a median metallicity (log(0/H)+12 = 8.30) that is 
~ 0.15 dex lower than from the targeted SN searches. 
This median metallicity is significantly different^ from 
the distribution of cither SN Ic (~ 0.3 dex lower median, 
Pks = l-7±f;° x 10~ 3 ) or the combination of the N = 54 

SN lb, Ic, and Ib/c (p K s = 3.1±l 2 7 9 x 10~ 3 ). As the 
SN Ic-BL metallicities derive almost exclusively from ex- 
plosion site spectroscopy, this difference is not dependent 
on spectroscopic methodology (pks = 7-2^g 2 Q 9 x 10~ 3 for 
SN Ic vs SN Ic-BL with explosion-site spectroscopy only) . 

It has been suggested that the metallicit y distribution 
of SN e Ic-BL may be bimodal (see e.g. IrVlodiaz et al.l 
l2008f ). This question is entangled with SN search 
methodology, because some SN Ic-BL are discovered via 
associated GRBs (by untargeted gamma-ray searches), 
while many SN Ic-BL without GRBs have been found by 
targeted SN searches. While we have shown SN Ic-BL 
are preferentially found at lower metallicity than other 
SN Ibc, Figure [7] illustrates that the metallicity distribu- 
tion of SN Ic-BL is broad, extending to super-solar metal- 

14 Statistical significance is verified despite a smaller sample size 
than in the SN lb versus Ic comparison because the difference in 
the distributions is much larger. Simulations of the type illustrated 
by Figurc[9]indicate that a sample of ~ 30 SN Ic+Ic-BL is sufficient 
to expose a 0.3 dex median difference at the pxs < 0-05 level in 
> 84% of trials. 
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licities. The Figure is visually suggestive of a bimodality, 
even among objects only from untargcted SN searches 
and with explosion-site metallicity measurements, but 
given the sample size the distribution is not significantly 
different than a Gaussian. Further studies comparing 
SN Ic-BL explosion properties (optical luminosity, pho- 
tosphcric velocity, etc.) with explosion site metallicity 
may inform the discussion of progenitor (sub) classes for 
these objects. 

4.4. SN lib vs. Ibc metallicities 

lArcavi et al.l (|2010l ) reported that SN lib preferentially 
occur in low-luminosity, likely low-metallicity, host en- 
vironments. However, when no cut is placed on SN 
discovery or spectroscopy characteristics, the difference 
between the combined sample distributions for SN lb 
(N = 47) and lib (N = 25) is not significant, with 
Pks = 0. 561^35 (see Figure [8]). Among only objects 
with explosion site spectroscopy (A^ib = 5 SN lib), 
the SN lib median log(0/H)+12 = 8.55, - 0.1 dex 
lower than SN lb, but this difference is not significant 
{pks = 0-32^0 24). Among SNc discovered by untargeted 
searches (iVnb — 8, mostly from this work), there is no 
suggestion of a difference between the host environment 
metallicity of SN lib and SN lb, with p KS = 0.75±g;||. 
As with SN Ic-BL, Figure [7] points to a bimodality in the 
SN lib metallicity distribution that cannot be statisti- 
cally verified with the present sample size. 

4.5. Type Ic supernovae from dwarf host galaxies 

Investigati ng SNe discovered in the first year of opera- 
tion of PTF, lArcavi et all 0fi$) found that SN lib and 
Ic-BL occur preferentially in dwarf host galaxies (defined 
by M r > — 18 mag), while they found SN Ic only in giant 
host galaxies (M r < — 18 mag). They interpret this as 
evidence for a dependence of the SN explosion properties 
on the metallicity of the explosion s ite. T ransforming the 
L - Z relation of iTremonti et al.l (|2004[ ) to the r band, 
they find that this threshold corresponds to a charac- 
teristic metallicity of 0.35 Z®, or log(0/H)+12 = 8.23 
on the PP04 scale. However, this photome tric approach 
has s everal limitations (see Section 3.1 of lArcavi et al.l 
120101) . In particular, the statistics are limiting, with 
only 6 SN Ibc in dwarf galaxies, and the absolute mag- 
nitude of the host galaxy does not necessarily reflect the 
metallicity of the explosion site accurately (Section l5.2p . 
Moreover, while the mass-metallic ity relation for star- 
forming galaxies is relatively tight ([Tremonti et al.ir2004l 
find a = 0.10 dex), the scatter is larger when lumino sity 
is used as a proxy for mass (ITremonti et al.l 120041 find 
a = 0.16 dex for Mb) and the scatter increases by a 
factor of > 2 at low- luminosities /mas s es in the dwarf 
regim e (see e.g. iKewlev fc Ellison! 120081 : iMannucci et al.1 
20111). Addit i onal scatter is introduced in the survey of 
Arcavi et al.l (|2010[) because they use a statistical trans- 
formation to convert their r-band magnitudes to the B- 
band and because they do not correct fo r extin ction 
of the host galaxies, which ITremonti et al.l (|2004f ) find 
is typically As ~ 0.3 mag ( ~ 0.05 dex i n Z) Fi- 
nally, the statistical analysis of lArcavi et al.l (|2010l ) de- 
pends on the arbitrary choice of the luminosity threshold 
(M r = -18 mag). 

Using the spectroscopic metallicity threshold 



log(0/H)+12 < 8.23, 4% of the combined sample 
of SN Ibc host galaxy observations we consider here 
are in dwarf hosts (6/171 with PP04 metallicity mea- 
surements). Of these six, two are SN Ic-BL discovered 
by untargeted surveys (2007ce and 2008iu, both with 
explosion-site metallicity measurements presented by 
this work), one is an untargeted SN lb (2007az; this 
work), one is a targeted SN lib (2 008ax; nuclear metal- 
licity from iKellv fe Kirshneill2011l) . and two are SN Ic. 
The SN Ic are SN2002jz, discovered by a targeted 
survey with an explo sion site metallicity measured by 
lAnderson et al.l (|2010D . and SN 2010Q, discovered by the 
untargeted CRTS survey and with explosion site metal- 
licity measured as a member of our Gold spectroscopic 
sample. However, the spec tral classificat i on of S N 2002jz 
is somewhat uncertain, as IHamuv et al.l (|2002bf ) suggest 
that it possibly displays Ha absorption reflective of 
Type lib SNe. Additionally, there are 3 SN Ic (all 
Gold classifications) in our sample for which we can 
not measure metallicities using strong line diagnostics, 
but are likely to be in sub-solar metallicity hosts and 
at least one of these (SN 20101z) is in an exceptionally 
low-luminosity host ga laxy (see Sect i on [51) . Finally, we 
note SN 2005kf, which iModiaz et all (|2011l ) report as a 
Type Ic SN from a Mb = —17 mag dwarf host galaxy 
with super-solar metallicity (log(0/H)+12 8.8). 

In summary, spectroscopic metallicity measurements of 
the host e nvironments of SN Ibc do not support the con- 
clusion of lArcavi et al.1 (|2010T l that SN Ic do not occur in 
low-metallicity host galaxies. This discrepancy may be 
due t o the small number of objects included in the sam- 
ple of lArcavi et al.l (|2010[) . Moreover, because secondary 
metallicity estimates made using host galaxy photometry 
introduce additional scatter, their role in detecting the 
subtle difference that could exist between the metallicity 
distributions of SN lb and Ic may be limited. In Sec- 
tion 22] we show that existing observations demonstrate 
that the difference in the median metallicity of SN lb and 
Ic host galaxies is almost certainly < 0.2 dex, while the 
systematic uncertainty in photometric metallicity mea- 
surements is ~ 0.2 dex. 

5. SYSTEMATIC EFFECTS 
5.1. Targeted vs. Untargeted SN Searches 

Galaxy-targeted supernova searches can bias SN host 
environment studies towards higher metallicities due to 
the galaxy L — Z relation, as illustrated in FigureQIJl The 
median metallicity measured for SNc discovered by tar- 
geted searches is log(0/H) +12 = 8.64, while for untar- 
geted searches it is log(0/H)+12 = 8.52 (-24% lower). 
The difference between the distributions is statistically 
significant (pks — 3 x 10~ 4 ). The fraction of galax- 
ies with metallicity < 1/3 Z© is larger by a factor of 
Njj/Nt ~ 5.1 in untargeted surveys, meaning that low- 
metallicity galaxies are strongly underrepresented in tar- 
geted surveys. This ratio is still appreciable for 1/2 Z & 
(Njj/Nt ~ 2.4), and even galaxies of solar metallic- 
ity are somewhat under-represented in targeted searches 
(Njj/Nt « 1-2). 

Consequently, galaxy-targeted SN searches offer a 
smaller baseline over which to probe for differences in 
metallicity distributions. The la spread in the metal- 
licity distribution of SN Ibc discovered by targeted SN 
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TABLE 10 

Statistics of the combined SN Ibc sample 





All Nuclear 


Exp. Site 




SN lb vs Ic 


All 

Targeted 
Untargeted 


47; 68; 0.08 j^gg 11; 26; 0.4812;^ 
26; 39; Omto,lt 7; 18; 0.57±°' |^ 
21; 29; 2.0^'* X 10" 2 4; 8; 0.38 f s 


36;42;0.16lo:?3 
19;21;0.63l°'38 
17; 21; 0.06^ 


SN lib vs lb 


All 

Targeted 
Untargeted 


25;47;0.56± ';35 13; 11; 0. 60±°; 3 g 
17; 26; 0.21^;^ 12; 7; 0.15±p " 
8;21;0.75l°;32 1;4 


12;36;0.3llE; : 4 
5; 19; 0.3210; 42 
7; 17;0.70lo; 2 5 


SN Ic-BL vs Ic 


All 

Targeted 
Untargeted 


18;68;1.3ti j x 10" 4 3; 26 
6; 39; 4.4±3 6 g 1 X 10" 2 2; 18 
12; 29; 1.7t?;° X 10~ 3 1;8 


15; 42; 1.8±;[ 2 5 3 X 10" 3 
4;21;0.38l°; 42 
11; 21; 7.2±J a 5 9 x 10~ 3 


SN Ic-BL vs Ibc 


All 

Targeted 
Untargeted 


18; 127; 2.3± 2 U { x 10" 4 3; 41 
6; 73; 4.7±4 3 5 1 x 10 ~ 2 2 ' 29 
12;54;3.ll2 2 7 9 x 10" 3 1; 12 


15;86;4.5l3 9 g 5 x 10" 3 
4; 44'; 0.39l°; 3 ? 
11; 42; l.ll^g X 10~ 2 



Note. — This table lists the KS test p- value (pks) an d sample size (AT) for SN Ibc 
mctallicity distributions (PP04 scale) in the combined sample. Each table entry reflects 
a different cut on the SN discovery (galaxy-targeted or untargeted) and spectroscopic 
methods (nuclear or explosion site; see also Figurc[8j . Each entry is given as JVi : N2 ; Pks ■ 
"SN Ibc" is a combination of SN lb, Ic, and Ib/c. 



searches in the combined sample is only 70% as large as 
that of untargeted surveys. As a result, the metallicity 
distributions of SNe are compressed and differences are 
reduced, as illustrated by Figure |U 

This L — Z bias has important consequences for studies 
which combine observations of SNe discovered by both 
targeted and untargeted surveys. To illustrate this, we 
take the assumption that SN lb and Ic share the same 
metallicity distribution and randomly draw samples from 
the targeted and untargeted SN Ibc metallicity distribu- 
tions shown in Figure HTJ1 We prod uce simulated samples 
constructed identically to that of iModiaz et al.l (|2011[ ) 
with respect to SN discovery characteristics: [6,11] tar- 
geted and [6,4] untargeted SN [Ib,Ic]. We find the SN Ic 
in the simulated samples to have higher average metal- 
licities than the SN lb in 64.2% of simulated trials, due 
solely to the L — Z bias. The average dif ference in the 
SN Ic and SN lb mctallicities reported by IModiaz et al.l 
(|201lD (> 0.14 dcx, depending on the diagnostic scale) is 
reproduced or exceeded in 3.3% of the trials. Similarly, 
we find that in 5.0% of trials, a KS test on the simu- 
lated observations would indicate that SN Ic host envi- 
ronments are significantly different (pks < 0.05) from lb 
environments. This represents a systematic effect biasing 
the results of the study, above and beyond the statistical 
ambiguity indicated by the KS test p- value. We conclude 
that differences in sample construction alone (the ratio 
of supernovae from targeted versus untargeted searches) 
can lead to erroneous differences in the metallicity distri- 
bution measured for the host environments of different 
SN types. 
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Fig. 10. — Cumulative mctallicity distribution for the host galax- 
ies SN Ibc (of all subtypes) discovered by targeted and untargeted 
searches, from the combined dataset (Section SNe with and 
without explosion site spectroscopy are shown (adopting only the 
explosion site measurement when both are available) ; if we restrict 
the sample to objects with explosion site spectroscopy, the effect 
is to shift to lower mctallicities by <C 0.05 dex. 

5.2. Isolating the SN Explosion Site 

Relative to nearby galaxies observed in targeted sur- 
veys, host galaxies in our exclusivcly-untargctcd sample 
will typically have smaller intrinsic radii and luminosi- 
ties and are typically found at greater distances. As a 
result, we are not able to resolve and measure the local 
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mctallicity at the explosion site of every host galaxy in 
our sample. As we discuss in Section 12.21 many of our 
measurements reflect integrated galaxy flux or only light 
from the brightest region of the galaxy (the nucleus). 
The situation is similar for SNe hosted by edge-on galax- 
ies, where the host environment spectrum is necessarily 
integrated over the full line-of-sight. 

It has been shown that SN Ic preferentially occur in the 
innermost regions of galaxies, m oreso than SN lb or lib 
(IPrieto et al.ll2008HKelly et al.| [200l lAnderson fc Jamesl 
20091; IHaberghamet al.ll2012D . In cases where we resort 
to measuring the galaxy nuclear metallicity, this implies 
that we would overestimate the metallicities of SN lb 
relative to Ic in the presence of a strong gradient (e.g. 
iZaritskv et al.|[T994D . We would therefore be less sensi- 
tive to the scenario where SN lb come from lower metal- 
licity environments than SN Ic within galaxies of similar 
nuclear characteristics. 

However, the role of gradients is not clear cut and the 
ability to resolve explosion environments is limited, even 
for surveys of local galaxies. Individual SN Ibc host envi- 
ronments (HII regions) can only be resolved from nearby 
emission nebulae for the most local events (e.g. z < 0.005 
for a 100 pc HII region and a 1" slit). Moreover, long- 
slit spectroscopy of nearby spiral galaxies illustrates that 
the oxygen abundance difference between the galaxy nu- 
cleus and disk outskirts is mo dest (very rarely as large 
as 0.3 dex; IMoran et al.l 120121 ). Metallicity gradients in 
galaxies should play an even smaller role for explosion 
sites c lose to the galaxy nucleus, where SNe Ibc usually 
occur (jKellv et al.ll2008l ). Moreover, nearby spiral galax- 
ies show sig nificant intrinsic scat t er in metallicity gra- 
dients (e.g. iFabbiano et al.l 120041 iRosolows kv fc Simon! 
120081 ISanders et al.M2012t but see also lBresolinll2011l ). 

We investigate the effect of physical resolution on the 
median metallicity difference measured between SN lb 
and Ic host environments in Figure [TT1 Monte Carlo sim- 
ulations incorporating the uncertainties in the individual 
metallicity measurements demonstrate that the differ- 
ence in the median metallicity of SN lb and Ic is never 
significantly larger than the uncertainty in the median 
for any sample. For comparison, we show a model where 
the true median difference of 0.2 dex that can be mea- 
sured accurately at z = 0, but is diminished linearly with 
increasing physical resolution until it is apparently dex 
at a modest rcdshift (z « 0.1, resolution of ~ 2 kpc). 
The model is shown cumulatively for a sample with the 
rcdshift distribution of the combined sample. 

Comparing the model to the existing datasets in Fig- 
urc[TT]suggcsts that physical resolution is not the limiting 
factor obscuring a difference in the median metallicity of 
SN lb and Ic. The median difference in the combined 
sample is not greater than 0.1 dex at any redshift (and 
always consistent with dex at 1 a), even though the 
corresponding model suggests that a 0.15 — 0.2 dex metal- 
licity difference would be preserved in the sample out to 
z > 0.2. We conclude that limitations due to physical 
resolution are not sufficient to mask a significant differ- 
ence between the metallicity distributions of the SN lb 
and Ic in the combined sample. 

A related issue is the methodological difference be- 
tween explosion site and galaxy-nucleus spectroscopy. In 
the combined sample (Section|4|), there are 30 SNe which 
have metallicity estimates at both the galaxy nucleus and 



explosion site. Figure [T2l indicates that nuclear metallic- 
ity measurements introduce both a small additional scat- 
ter (rms discrepancy of 0.13 dex) and bias (median differ- 
ence of —0.08 dex) relative to explosion site spectroscopy. 
The sign of the effect is consistent with offset explosion 
sites being systematically lower-metallicity environments 
than galaxy nuclei due to galactic metallicity gradients, 
with the effect strongest among the highest-metallicity 
host galaxies (nuclear Z > Zq). In the median the dis- 
crepancy between the nuclear and explosion site metallic- 
ity measurements is only 0.68 times the uncertainties in 
the individual metallicity measurements. We do not find 
evidence for a significant difference between the explosion 
site versus nuclear metallicity discrepancies of SN lb and 
Ic in the combined sample; their medians arc different by 
only 0.03 dex, which is much smaller than the rms. 

Simulations of the type illustrated in Figure [9] indicate 
that the bias and scatter introduced by nuclear spec- 
troscopy should be small for studies of SN Ibc host en- 
vironments. Introducing an additional 0.1 dex scatter 
for all measurements, the sample size necessary to re- 
cover a 0.1 dex metallicity difference between SN lb and 
Ic increases by only a small amount (~ 20%). Intro- 
ducing a 0.1 dex bias only for objects with super-solar 
metallicity (and assuming intrinsic median metallicities 
of log(0/H)+12= 8.5 and 8.6 for SN lb and Ic, respec- 
tively) , we find a negligible effect on the significance indi- 
cated by the KS test. In summary, the combined sample 
demonstrates that nuclear metallicity measurements can 
be used as effectively unbiased tracers of SN Ibc progen- 
itor metallicity. 

A third issue relating to the isolation of the SN explo- 
sion site is the use of ph otometry to estima te the nuclear 
galaxy metallicity (e.g. lArcavi et al1l2010l ). or addition- 
ally using a simple metallicity gradie nt model to estimate 
the e xplosion site metallicity (e.g. iBoissier fe Prantzosl 
120091) . In either case, photometric mctallicity estimates 
will carry the same scatter as spectroscopic nuclear 
metallicity measurements (~ 0.1 dex, see above) and an 
additional uncertainty > 0.16 dex due to the scatter in 
the L — Z relation (see Section |475|) . Assuming a repre- 
sentative 0.2 dex uncertainty for photometric mctallicity 
estimates and performing the same simulations as in Fig- 
ure^ we find that a small (0.1 dex median) metallicity 
difference between SN lb and Ic host environments can 
only be detected with a sample ~ 2 times as large as 
an equivalent sample of spectroscopic metallicity mea- 
surements. Our conclusion that the combined sample of 
> 100 SN Ibc has not revealed a difference between the 
mctallicity distributions of SN lb and Ic m a y the refore 
conflict with the suggestion by lArcavi et al.l (|2010f ) that 
a sample of 18 photometric metallicity measurements for 
SN Ibc host galaxies could reveal a significant difference. 

In summary, the role of explosion site spectroscopy in 
SN host environment studies is complex. The ability to 
isolate SN explosion sites is limited for smaller galaxies 
and at larger redshifts (e.g. those found by untargetcd 
surveys). However, given the significantly larger base- 
line for making measurements of metallicity distributions 
offered by untargeted SN searches (Section 15. ip . the ef- 
fect of physical resolution is likely to play a secondary 
role. Substituting nuclear for explosion site spectroscopy 
can lead to discrepancies in metallicity of ~ 0.1 dex, 
but we find that small differences in mctallicity distri- 
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Fig. 11. — Bottom: the cumulative difference measured between the median mctallicity of SN lb and Ic host environments at rcdshift < z 
where > 4 SN lb and Ic are available, using log(0/H) as estimated from the PP04 diagnostic. The black dashed line illustrates the 
model for the effect of physical resolution on the measurement of a 0.2 dex difference in the me dian , described in the text. Except for the 
K+11 sample, only objects with explosion site spectroscopy claimed are considered (see Scction l4.ll) , The shaded region illustrates the la 
bootstrap errors on the median difference in the combined sample. Top: redshift distribution of SN lb and Ic in the individual studies of 
the combined sample (see Table Q5). 

buttons could still be recovered by studies employing nu- 
clear spectroscopy. However, the additional uncertainty 
introduced by photometric mctallicity measurements is 
larger and may limit the role of these investigations in 
discovering differences in metallicity distributions, unless 
sufficient statistical power can be gained by larger sample 
size. 



5.3. Uncertainties in SN classification 

The spectral classification of SNe could also potentially 
effect the metallicity distribution inferred for the host 
environments of each SN type (see also Section l2~Tj) . 

It is unlikely that a SN Ic would be mistaken for a 
SN lb, if the type is designated based on a clear detection 
of He I lines, but it is possible for a SN lb to be classified 
as an SN Ic if the He lines are weak and S/N is poor, or if 
the He lines have not yet developed at the epoch of spec- 



troscopy (jHamuv et al.ll2002a| ). Misidentifying SN lb as 
SN Ic could act to blur any distinctions in the metallicity 
distribution of their host environments. We evaluate the 
possibility that this blurring could hide a large (0.2 dex) 
difference in the metallicity distributions of SN lb and Ic 
host environments for a dataset as large as the combined 
sample using simulations. We randomly draw SN lb and 
Ic metallicities from Gaussian distributions with widths 
of 0.2 dex and then mis-attribute a certain fraction of the 
draws to SN Ic. We find that the KS test indicates a sta- 
tistically significant difference (pks < 0.05 at la) unless 
a large percentage (> 2/3) of SN lb are mis-classified as 
SN Ic. 

It can be difficult to distinguish normal SN lb from 
SN lib, especially if spectra are not available near peak. 
There is a diversity in the hydrogen mass surrounding 
SNe lib ([Chevalier fc Soderberd[2010() and, spectroscop- 
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Fig. 12. — Comparison of PP04 metallicity measurements for 
30 SN Ibc host environments from galaxy nucleus and explosion 
site spectroscopy, from the combined dataset. The horizontal line 
marks equality and the dashed vertical line marks the solar metal- 
licity. The colors reflect SN type as in Figurc|8]and the sizes encode 
the rcdshift of the object. 

ically, SN lib seem to r epresent a transitional t ype be- 
tween _SNJbaj^NIlj|Mni^ How- 
ever, lArcavi et al.l ([2010D have suggested that SN lib 
prefer low-mass, low-metallicity host galaxies. It is dif- 
ficult to reconcile these observations into a consistent 
progen itor model, but the conclusions of lArcavi et al.l 
(I2010D imply that gr oupi ng SN lb and lib to gether (as 
iLeloudas et al.l 120111 and iModiaz et al.l 120111 have done 
explicitly) could bias downward the metallicity distribu- 
tion inferred for SN lb environments. 

Similarly, if SN Ic-BL are grouped together with SN Ic, 
the SN Ic metallicity distribution could also be biased 
downward ( Section 14. 3|) . There is no clear definition of 
the ejecta velocity or line width that distinguishes SN Ic- 
BL from SN Ic, and any distinction will strongly depend 
on the spectroscopy epoc h due to the veloc i ty evolution of 
the p hotosphere (see e.g. IPian et ail feOOB; Sa nders et al.l 

[ml. 

We have attempted to reduce the effects of SN 
classification errors by revisiting the SN spectroscopy 
on a case-by-case basis, rather than relying on the 
classification reported in the IAUCs, and by separat- 
ing our dataset based on our degree of confidence 
in the classifications (Gold and Silver samples; Sec- 
tion 12.11) . Of t he previous stud i es in our c ombined sam- 
ple ( S ection |4|, ILeloudas et al.l (|2011[ ) and IModiaz et al.l 
(120111) made similar re visi ons to SN classifications, while 
lAnderson etHI ([2010f ) and lKellv k Kirshnerl (pOllh have 
relied on classifications reported in the circulars. 

5.4. Selection effects in spectroscopic follow-up 

The number of optical transients detected by SN 
searches often exceeds the resources available for spec- 
troscopic follow-up. For example, because SNe with- 
out host galaxies clearly detected in discovery images 
may be mistaken for non-SN optical transients, the 
spectroscopic follow-up for some SN searches will be 
biased against SNe with low-luminosity host galax- 



ies (|Drake et al.l l2009al ) . Similarly, some untargeted 
SN searches such as ESSENCE and SDSS-II have fo- 
cused on disc overing SNe la at the exclusion of core- 
collapse SNe (|Sako et al.|[200l iFolev eTail [2009) . Be- 
cause the host galaxy properties of SN l a differ sys- 
tema t ically from those of S N Ibc (see e.g. IKellv et all 
120081: iMannucci et al.l 12001) . this selection effect could 
influence resulting studies of SN Ibc host environments. 
Moreover, some optical surveys are less sensitive to SN 
detection in the central regions of galaxies. However, 
because the photometric and host galaxy properties of 
SN lb and Ic arc significantly more similar to other sub- 
classes of SNe, the systematic effects of biases in SN 
search spectroscopic follow-up is likely to effect SN lb 
and Ic in similar ways. Therefore we expect it to be a 
second order effect influencing the results of this study, 
likely much smaller than the previous effects discussed. 

In a magnitude limited survey, th e ratio of SN Ic to 
lb has been estimated to be ~ 1.6 (|Li et al.l 120 111) . In 
our study, this ratio is similar, 1.8, suggesting no sig- 
nificant bias for or agai nst SN lb or Ic . For compari- 
son, this ratio is 1 .4 for lAnderson et al.l (I2010D. 2.1 fo r 
IKellv k Kirshnerl (l20ll. 5 for ILeloudas et al.l (I201TI) . 
and 1.25 for IModiaz et all (pOlTT ). 

SN Ic-BL, comprising 13% of our sample, are over- 
represented in our sample with respect to galaxy- 
targeted, volume-limited surve ys (e.g. < 5% of SN Ibc 
in LOSS, not including SN lib. ILi et al.|[2^Tl . This is a 
natural consequence of two factors. First, SN Ic-BL typ- 
ically have brighter pe ak luminosities tha n other SNe Ibc 
(~ 1 mag in i?-band. iDrout et al.ll2011f ) and can there- 
fore be discovered over a larger volume in a magnitude- 
limited survey. Indeed, the median redshift of SNe Ic- 
BL in our sample (z = 0.056) is nearly twice that of 
the other SNe Ibc in our sample (z = 0.034). Second, 
we have found that Ic-BL preferentially occur in low- 
metallicity galaxies that will have preferentially lower lu- 
mino sities (see also lArcavi et al.ll2010T : IKellv k Kirshnerl 
1 2 1 If ) . These galaxies are therefore under- represented in 
galaxy-targeted SN searches, but are not excluded from 
the untargeted SN searches we draw our sample from. 
We further note that the percentage of SN l ib in our sam- 
ple (1 7%) is similar to the value found bv lSmartt et al.l 
((2009T) among those SN classified and reported in the I AU 
circulars: ~ 16% of SN Ibc. However. ILi et al.l (|201lD re- 
port an SN lib rate > 2 t imes this value based on SNe 
found by LOSS. ILi et al.l (pOll that the LOSS classifi- 
cations are based primarily on photometry (in particu- 
lar, identification of the double-peaked SN lib lightcurve 
shape), while classifications reported in the circulars are 
primarily based on single-epoch spectroscopic observa- 
tions. Finally, we note that any effect due to the cosmic 
star formation history should be negligible over the mod- 
est re dshift range (z < 0.3) of our sample (jGrieco et alJ 
120121) . 

An additional systematic effect could act if the host 
galaxy properties of SN Ibc are found to correlate with 
the exp losion peak magnit ude, as has been found for 
SN la (jHamuv et all fl996f). For example, if it were 
the case that SN Ic-BL in low-mass/metallicity galaxies 
are brighter than their counterparts in brighter galaxies, 
they would be over-represented in samples of host galaxy 
spectroscopy. This potential effect could be evaluated 
more thoroughly with more extensive studies of SN Ibc 
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lightcurve pro perties, as the sample size of existing stud- 
ies are small ([Richardson et al.ll20"ol iDrout et all 120111: 
iLi et al.1l2011l) , or by performing a volume-limited survey. 
Moreover, dust obscuration co uld prevent the disc overy 
of a significant fraction of SNe (jMattila et al.l[2012l ) , and 
could potentially be correlated with both host galaxy and 
explosion properties. 

5.5. Depth limits for host galaxy spectroscopy 

Due to the galaxy L — Z relation, the lowest metallicity 
galaxies will also have the lowest luminosities, and may 
therefore be under-represented in spectroscopic studies 
that require sufficient S/N in the nebular emission lines 
to derive metallicity. Table [6] indicates that we are 
potentially excluding 3 SN Ic in low-mctallicity hosts 
from our sample due to insufficient S/N, while we are 
likely only excluding one such SN lb. The effect of 
these exclusions is significant: if we suppose that each 
of the excluded galaxies falls at the low-mctallicity end 
of the observed distribution, than the median mctallic- 
ities we measure for the SN [Ib,Ic] in our sample be- 
come log(0/H)+12= [8.46,8.42]. Considering these ex- 
cluded galaxies indicates that the median metallicity for 
SN Ic host environments could be lower than for SN lb 
in our sample. Such a revision would not conflict with 
the results of our statistical tests, which indicate that 
the difference we measure in the metallicity distribu- 
tions of SN lb and Ic host environments is not significant 
(Section 13. 4[) . This illustrates the role of small number 
statistics in spectroscopic studies of SN Ibc host environ- 
ments. Moreover, it advocates for continued follow up of 
supernova host environments with facilities capable of 
measuring metallicity for the low-luminosity, relatively 
high-redshift host galaxies discovered by untargeted SN 
searches. 

6. DISCUSSION 
6.1. SN lb and Ic progenitor models 

The new observations presented in this work, and pre- 
vious observations synthesized in our combined dataset, 
suggest and constrain differences in the progenitor star 
population of SNe Ibc. The modest difference observed 
in the median metallicity of SN lb and Ic host environ- 
ments corresponds to a very small difference in mass- 
loss rates for single star progenitors. In general, the 
mass loss rates of massive star s is taken to scale as 
M ~ {Z/Z Q ) m where m = 0.5 dKudritzki et all [l98l 
or pe rhaps m = 0.86 for WR winds (jVink fe de Koten 
l2005|. The difference in median metallicity we measure 
for SN lb and Ic from explosion site spectroscopy of Gold- 
classification SNe in our untargctcd study (Section l3.3j) 
then corresponds to a difference in mass loss rates of a 
factor of Mjc/Mn, ~ 1.17 (1-32) for the power law slope 
m = 0.5 (0.86). If we use instead the median difference 
inferred from the combined sample using only explosion- 
site spectroscopy ( Section |4~2|) . the mass loss rate differs 
by a factor of Mi c /Mib ~ 1.11 (1.20). In this simplistic 
analysis, it seems improbable that this small difference 
in mass loss rate is sufficient to strip the entire He layer 
from the progenitor star in order to produce the envelope 
composition indicated by the spectrum of the explosion. 

In more detail, the difference in the SN lb and Ic 
metallicity distributions can be interpreted in terms of 



a metallicity-dependence for the critical initial-mass re- 
quired for a progeni t or sta r to explode as SN lb or Ic. 
iBoissier fc Prantzosl (|2009t) estimated how this critical 
mass depends on metallicity by comparing the observed 
difference in the rates of SN lb and Ic in different metal- 
licity bins to a model that has an explicit dependence 
of the threshold mass for SN lb and Ic explosions. Us- 
ing photometric metallicity estimates, they find that the 
critical mass varies by a factor of ~ 2 over a factor of 
~ 3 in metallicity. However, our combined dataset in- 
dicates that a difference in the metallicity distribution 
of SN lb and Ic has yet to be measured robustly (Sec- 
tion S2]), suggesting that the metallicity-dependence of 
this critical mass may be much more subtle. 

Moreover, the observed relation between galaxy mass, 
metallicity, and star formation rate implies that metal- 
poor galaxies typically have higher specific star forma tion 
rates (jLara-Lopez et al.ll20T 0: Mann ucci et al.l[2010( ). In 
a single star progenitor model where SN Ic are produced 
by more massive stars, this could indicate an elevated 
rate of SN Ic relative to SN lb in metal-poor galax- 
ies if the star formation events are short-lived. This 
could potentially mask the effect of metal-line depen- 
dent winds. Our measurements of the ages of young 
stellar populations of SN Ibc host galaxies do not in- 
dicate that SN Ic come from younger stellar populations 
(Section I3.5[) . but a larger sample is needed to address 
this question in detail. A metallicity-dependent slope for 
the initial mass function c ould also effect SN Ibc rates 
(jPrantzos fe Bo issicr 20031) . Assuming SN Ic progenitors 
are more massive than those of SN lb, a top-heavy IMF in 
low- metallicity environments would elevate the SN Ic/Ib 
rate and could push the metallicity distribution of SN Ic 
d ownward. 

iSmartt et al.l <|2009f ) suggest that luminosity limits 
from pre-explosion imaging and ejecta masses from light- 
curve modeling rule out massive WR progenitors for a 
sample of nea r by SN Ibc. However, the SN Ibc studied by 
ISmartt et al.l (I2009T ) are primarily in nearby, high-mass 
host galaxies that are likely to reflect high-mctallicity 
progenitor stars. In essentially all progenitor models, 
SN Ibc progenitors are expected to be less massive (and 
less luminous) at higher metallicities. Indeed, the rate 
of SN Ibc in nearby (z < 0.04), low-metallicity envi- 
ronments is expected to be quite low (at best a few per 
year) , makin g direct progenitor detection observationally 
challenging ( Young et al.l 1200*81 ) . Moreover, local extinc- 
tion may play a larger role in o bscuring SN progenitor 
stars than previously recognized (|Walmswell fc Eldrided 
[2011). 

iPodsiadlowski et al.l (| 1992T) suggested that Roche lobe 
overflow via binary interaction could be responsible for 
stripping significant amounts of material from the pro- 
genitor stars of SN Ibc. Observations of OB stars indi- 
cate that a significant percentage of potential Type Ibc 
SN prog enitors are likely to be in interacting binary 
systems (iKobulnickv fc Frverl 120071: iKouwenhoven et al.l 
120071 ISana et al.l 12012^ lEldridge et al.l (|2008fl have 
shown that the observed metallicity-dependence of the 
relative rates of SN II and Ibc can be reproduced using 
binary population synthesis models informed by the ob- 
served populations of red supergia n ts, W olf-Rayet, and 
other massive stars. iSmith et al.l ([201 ll ) have argued 
based on SN rates that the majority of SN Ibc pro- 
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genitors may come from binary star systems. However, 
additional modeling is needed to predict the relative 
metallicity-distribution of SN lb and Ic that would re- 
sult from binary progenitor stars. Moreover, the effect 
of mixing may further complicate the comparison of ob- 
servations to SN Ibc progenitor mod els for both single 
and binary stars (|Dessart et al .1120121) . 

6.2. Comparison to nearby GRB-SNe 

In the "collapsar" model, the progenitors of LGRBs 
found in association with SNIc-BL (GRB-SNe) arc 
massive stars with high rates of core rotation, imply- 
ing sub-solar metallicities [Z < 0.3 Z@) in order to 
minimi ze angular momentum l osses due to line-driven 
winds dWooslev fc Bloomll2006f ). The observational re- 
sult that most LGRBs and GRB-SNe are discovered 
in dwarf, sub-solar metallicity galaxies has been inter- 
prete d as evidence suppor t ing this model ( Fruchter et al.l 
[200l iStanek et all [200l iLevesaue et all l2010al) . Be^ 
cause SNIc-BL have traditionally been found in higher- 
metallicity environments, a "cut-off" metallicity has 
been proposed to distinguish stars which will produce 
GRB-SNe from those that w i ll produce non-relativ istic 
SNIc-BL (|Modiaz et al.ll200l iKocevski et al]l200flh . 

However, several discoveries have challenged the role 
of metallicity in the production of LGRBs. LGRBs 
and rclativistic, engine-driven SNe have been found 
in super-solar metall i city en v ironments (iBerger et ail 
20071 ILevesaue et al.) I2010dt iSoderberg et al l I2010bt 
Levesque et al.l I2010bt iGraham et all 12009( 1. 



SNIc-BL with strong limits on the association of rela- 
tivistic ejecta have be en found in sub-sola r metallicity 
environments (see e.g. I Sanders et al.l [20 11(1 . This shift 
in the discovery environments of GRBs and SNIc-BL 
can be explained in terms of the transient search strat- 
egy — GRBs are discovered via their gamma-ray emis- 
sion through untargeted searches, while past studies of 
SNIbc host environments drew primarily from galaxy- 
targeted searches, which are heavily biased towards 
higher- metallicity host environments (see Section 15. ip . 
The advent of wide-field, untargeted SN searches has 
enabled the discovery of SNIc-BL in sub-solar metal- 
licity env^romnents^Indjed, the results from this work 
and iKellv fe Kirshnerl (|2011l ) demonstrate that SNIc- 
BL preferentially occur in lower-metallicity environments 
than SNIb or Ic (see Section 1331 and Section |4~3|) . 

Alternatively, the observed metallicity distribution of 
LGRB host environments has been interpreted as a 
secondary manifestation of a preference for high-star 
formation rate envir onments (|Kocevski fc Westl 120111 : 
iMannucci et al.ll201lH . Our results are consistent with 
this view. In Section 13751 we found that SN Ic-BL host en- 



vironments potentially have younger stellar populations 
than those of other SNIbc and more similar to nearby 
LGRBs. In Section l3~3l we find that the untargeted SN Ic- 
BL in our sample have a median metallicity ~ 0.2 dex 
higher than that of the nearby LGRBs (z < 0.3), but 
given the small size of both samples (N = 7 and 7, re- 
spectively) , the KS test indicates no significant difference 
between the distributions (pks = 0.13). 

In summary, observations of SNc discovered by untar- 
geted searches have substantially reduced the discrep- 
ancy reported between the metallicity distribution of 
SNIc-BL with and without associated LGRBs, while ver- 



ifying a shared preference for environments with high- 
star formation rates and/or very young stellar popula- 
tions. These findings are consistent with the view that 
massive stars (M > 40 Mq) are the common progeni- 
tor for both these types of explosions, but suggest that 
metallicity does not play the primary role in the for- 
mation of a central engine. Similarly, ILevesaue et al.l 
(|2010eH have found that host environment metallicity 
does not correlate to t he ga mma-ray energy release of 
LGRBs. iGeorev et all (1201211 ha ve suggested, using the 
models of lEkstrom et aLl i 2012( 1. that differential rota- 
tion in LGRB progenitors could moderate the coupling 
between WR winds and the stellar core and may be re- 
sponsible for reducing the role of metallicity in the ex- 
plosion. 

6.3. Studies of SN Ibc in the LSST era 

In the coming years, the discovery rate for SNe 
promises to grow dramatically as existing high-cadence, 
wide field surveys continue to operate and new surveys 
such as the Dark Energy Survey and the Large Synoptic 
Survey Telescope (LSST) come online. As LSST could 
discover ~ 10 6 SNe per year, it is anticipated that the 
SN discovery rate in this era will far outstrip the ca- 
pacity t o perform spectroscopic follow-up on ind ividual 
objects (|LSST Science Collaborations et all 12009(1 . This 
significant increase in the SN discovery rate could per- 
mit studies of untargeted SN Ibc host environments with 
sample sizes sufficient to detect even small (< 0.1 dex) 
differences in the SN lb and Ic metallicity distributions 
without statistical ambiguity. We suggest that a sample 
of > 100 such o bject s would be required for that pur- 
pose (see Section [4~2j) , but caution that attention must 
continue to be paid to the systematic effects discussed in 
Section [5] 

Moreover, such an SN discovery rate suggests that SN 
classification will become the limiting factor for certain 
science goals. Investigating the GRB-SN connection will 
be more difficult if a sample of well-classified SN Ic-BL 
cannot be identified. In SN la samples for precision 
cosmology that are assembled using light curve-based 
classification alone, SN Ibc may be major contaminants 
(~ 20 %) because their light curve shape resembles that of 



SNIa (IGong et aIJ201ttlKessler"eTaiTMol:lGiergo et al.l 
120121 ). The potential for ambiguity between SN la and 
SN Ibc light curves is especially great in the case of SN Ic- 
BL, because the average absolute peak magn itude is more 
simil ar to SN la than for normal SN lb or Ic (jDrout et al.l 

We therefore suggest that information about the dis- 
tribution of host galaxy properties for SN la and Ibc 
could be used to inform photometric SN classification in 
the LSST era. Because SN Ic-BL preferentially occur 
in lower luminosity, low er metallicity host galaxies than 
other SN Ibc or SN la (|Li et al.ll20lil ; IKellv fc Kirshnerl 
1201 11 1, the difference in flux between the SN and the host 
galaxy could be used to distinguish SN Ic-BL from other 
SN types for the dual purposes of 1) identifying SN Ic- 
BL for follow-up SN studies and 2) filtering SN Ibc from 
samples of SN la for precision cosmology. However, a 
study using information about the host galaxy distri- 
bution to inform SN classification would sacrifice the 
ability to make unbiased inferences about the proper- 
ties of the host galaxy populations (similar to the case of 
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galaxy-targeted SN surveys) and could introduce system- 
atic effects if used to construct samples for SN la cosmol- 
ogy. Moreover, it would also exclude events with unex- 
pected properties, such as the engine-driven Type Ic-BL 
SN 2009bb t hat exploded in a lumi nous, high-mctallicity 
host galaxy (jLevesque et al.lfeOlOd ). 

7. CONCLUSIONS 

We have presented the largest study to date of 
SN Ibc host environments unbiased with respect to the 
galaxy L — Z relation as imposed by galaxy-targeted SN 
searches. By reporting mctallicity measurements for 50 
objects, we more than double the number of host envi- 
ronment metallicity measurements for SN Ibc discovered 
by untargeted surveys. 

We conclude that: 

1. SN Ibc host environments discovered through 
targeted surveys are significantly biased towards 
higher mctallicitics, representing host environ- 
ments that are typically more enriched by ~24%. 
The ratio of low- mctallicity (Zq/3) host environ- 
ments probed by untargeted versus targeted SN 
searches is Njj/Nt ~ 5.1, and galaxy-targeted SN 
searches offer a smaller baseline for probing metal- 
licity distribution differences. 

2. In our own sample and combining observations 
from all spectroscopic studies of SN Ibc host envi- 
ronments to date, we find no statistical difference 
between the metallicity distributions of SN lb and 
Ic or between SN lb and lib. We place a limit on 
the median metallicity difference between SN lb 
and Ic (< 0.1 dex) and find that a sample > 2x as 
large would be required to unambiguously confirm 
a difference at that level. This limit corresponds 
to a very small difference in the mass loss rate of 
metal line-driven winds, suggesting that it may not 
be the dominant factor distinguishing SN lb and Ic 
progenitors. 

3. SN Ic-BL are found in host environments with sub- 
stantially low er metallicity than SN Ic , confirming 
the result of iKellv fc Kirshnerl (|2011[ ). We show 
that the median metallicity of SN Ic-BL found by 
untargeted SN searches is ~ 0.15 dex lower than 
those found by targeted SN searches, yielding closer 
agreement to LGRB host environments. Moreover, 
the young stellar populations of SN Ic-BL appear 
to be lower than SN lb and Ic, but the sample is 
not large enough to be significant. 

4. Evaluating the systematic effects afflicting stud- 
ies of SN Ibc host environment metallicities, we 
find that the bias introduced by galaxy-targeted 
SN searches is most significant. Galaxy- nucleus 
spectroscopy can serve as a good proxy for explo- 
sion site mctallicity measurements, but studies us- 
ing photometric mctallicity estimates would require 
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sample sizes ~ 2 times larger. Biases in spectro- 
scopic follow-up of SNe discovered by optical tran- 
sient searches and uncertainties in SN classification 
play smaller roles. 



In the era of wide-field, untargctcd SN searches, we 
anticipate that significant progress will be made towards 
unveiling the progenitors of SN Ibc through the study 
of their host environments. We advocate for continued 
spectroscopic studies of the host environments of SN Ibc 
discovered by untargeted surveys to uncover or place 
stricter limits on the difference in the metallicity distri- 
bution of SN lb and Ic. Optical facilities with large light- 
gathering power are required to measure metallicities in 
the low-luminosity, relatively high-redshift galaxies host- 
ing SNe discovered by untargeted surveys. To identify 
SN Ic-BL in this era and make progress on the GRB-SN 
connection in the face of limited resources for spectro- 
scopic follow up, we suggest that host galaxy properties 
will become a significant aid to future photometric SN 
classification schemes. 
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